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Abstract—The dispersed phase holdup for spray columns can be related to the phase flow 
rates by the following equation, derived originally for packed columns [5] and applied subse 
quently to pulsed [10] and rotary annular columns [11] : 
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In this expression, v, is the characteristic velocity identified as the mean relative velocity of 
the droplets extrapolated to essentially zero flow rates. 


Differentiation of the above equation, treating x as the independent variable, gives expressions 
which have been used to correlate flood point data in terms of v, [10, 11]. It has also been 
shown that the limiting holdup at the flood point is dependent only on the phase flow ratio and 
is independent of the physical properties of the system. 


In the case of spray columns, it has been found possible to predict vy on the basis of the normal 
correlation of the free falling velocity for spheres, assuming the droplet size to be given by 
Haywortn and TreyBav’s relationship [6]. The entire flooding curve can thus be predicted 
for any system in terms of its physical properties and the dispersed phase nozzle geometry. 


Résumé—Les auteurs ont montré que le hold up de la phase dispersée dans les colonnes a 
pulvérisation peut étre relié & son débit par l'équation suivante, obtenue originellement pour 
des tours 4 garnissage et appliquée par la suite aux colonnes pulsées et & disques tournant : 
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Dans cette expression, vg est la vitesse caractéristique identifiée a la vitesse moyenne relative 

des gouttelettes extrapolées a des débits rigoureusement nuls. 


La différentiation de léquation ci-dessus, 2 étant considéré comme variable indépendante, 
donne des expressions utilisées 4 relier le point d’engorgement en fonction de Up: Ils ont aussi 
montré que le hold up limite au point d’engorgement ne dépend que du rapport des phases et 
est indépendent des propriétés physiques du systéme. 


Dans le cas des colonnes a pulvérisation, les auteurs ont montré qu’il était possible de calculer 
Up, en se basant sur les vitesses normales de chute libre de sphéres, avec 'hypothése que la grosseur 
dex gouttes est donnée par la relation de Haywortn et TreyspaL. La courbe d’engorgement 
toute entiére peut alors étre calculée pour un systéme quelconque au moyen de ses propriétés 
physiques et de la forme géométrique de l'ajutage de dispersion. 
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1 INTRODUCTION 
IN rue design of spray towers for liquid-liquid 
extraction operations, a knowledge of the limiting 
flow rates of the two phases is necessary before 
the column diameter can be calculated, and with 
this object in view a number of equations have 
been proposed for correlating flooding data in 
Thus CoLtpurn [3] 
plotting V, Ap against Ve V as 
but this 


successful in bringing together the results obtained 


spray towers. suggested 
abscissa on 
logarithmic co-ordinates, was not 
Some improvement could 
be obtained by dividing the ordinate by the 


of the droplet diameter, but the 


for different systems. 
square root 
correlation was not good, 

A semi-theoretical equation, incorporating the 
droplet diameter and two empirical constants, 
was proposed by BLanpineG and Exner [1], but 
as pointed out by these authors, their relationship 
was not an exact representation of the flooding 
mechanism and required some labour in its 
manipulation, Exar and Foust [4], using solid 
spheres instead of deformable liquid droplets, 
obtained data which paralleled those of BLANDING 
and Every, and a method of correlation was pro- 
posed which was based upon a knowledge of the 
particle diameter and upon the assumptions that 
the same factors are operative in spray towers 
and fluidized beds, and also that the void fraction 

Reynolds number relation of Witnetm and 
Kwavk [13] is the same for the two. However, 
calculations of the holdup at the flooding point, 
based on the latter assumption, gave values in 
the neighbourhood of 0-05 [12], whereas experi- 
mental determinations are usually in the range 
0-10-0-45 and can, in theory, range from 0-0.50, 
as will be shown later. 

Mixarp and Jonnson [8], in an extensive 
study of spray column behaviour using a variety 
of aqueous-organic systems, obtained a semi- 
theoretical equation involving the droplet ‘size 
and the dispersed phase holdup at flooding, based 
Bertrertri’s theory for 

In order to translate 


upon an extension of 

gas-liquid packed towers. 
the equation into a more usuable form, the droplet 
size and the holdup were replaced by an experi- 
mentally determined “ nozzle constant.” How- 
ever, the necessity of having to determine this 


constant for each particular nozzle geometry 
imposes an immediate limitation on the method, 
Another correlation involving droplet size and 
holdup at flooding has been proposed by Sakiapis 


and Jounson [9] who, in an attempt to simplify 


1 O75 
on the 
22 


grounds that this quantity tends to approach 


their result, dropped the term | l 


unity at high holdups. However as pointed out 
previously, the fact that a only tends to 0.50 in 
the limit renders this procedure questionable. 
Of the methods available for the prediction of 
holdup at the flooding point, that proposed by 
Witnecoa and Kwack has been mentioned above, 
but was discounted on the grounds that it gave 


Both Minarp 
and JOuUNSON and SAKIADIs and JouNSON derived 


values which were on the low side. 


expressions for the holdup in terms of the limiting 
flow rates, the densities of the phases and the 
droplet size. However, the latter expression which 
is derived for true flooding conditions as opposed 
to total rejection indicates a fractional holdup of 
unity at zero continuous phase flow, which is 
contrary to experimental evidence [8], the limiting 
In this 
connection BLANDING and ELGIn reported that 


value being 0-50 as will be shown below. 


“in no case did the holdup tend to equal in 
volume that of the column, or to fillit. A tendency 
for the discontinuous phase to coalesce and be- 
come continuous, filling the column at flooding, 
was found only with end designs which produced 

the entrance to the 
seen that no 


marked coalescence at 


column.” It will be therefore, 
entirely satisfactory method of estimating the 
limiting holdup has yet been proposed. 

Since previous work had shown [1, 8] that spray 
columns could be operated up to the flooding point 
with little coalescence of the dispersed phase, it 
seemed reasonable to suppose that the data could 
be correlated by the same method as that used 
previously for pulsed and rotary annular columns 
[10, 11]. In the present paper, equations are 
presented which enable the limiting holdup to be 
calculated from a knowledge of the flow ratio only 
and, used in conjunction with Hayworrn and 
TREYBAL’s equation [6] for droplet size, can be 
used to predict the effect of variations in physical 
properties and nozzle diameter on the flooding rate. 
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2. THEorRETICAL TREATMENT 


(a) Limiting Flow and Holdup Equations 


It has been shown [5], that for packed columns, 
the dispersed phase holdup can be related to the 
phase flow rates by the equation : 

V, 


«VQ r) (1) 


1 a 

where I’, is the characteristic velocity, identified 
as the mean relative velocity of the droplets 

Plots of this 
equation in the form V, and V, lagainst x for a 
typical value of V, = 500 ft hr are shown in 
Fig. 1. Although equation (1) is cubic, only the 
root which lies between zero and the flooding 


extrapolated to zero flow rates. 


point values of Vand V. is realizable in practice. 
These portions of the holdup curves are indicated 
by the full lines in Fig. 1. 
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Fig. 1. Plot of Eq. (1) in the form Vy and JV, against 


for a typical value of Vy = 500 ft/hr. 


Equation (1) was found [10, 11] to be applicable 
up to the flooding point for mechanical columns, 
because in these cases coalescence is reduced to 
a minimum and the droplet size distribution for 
a given system is a function of the column geo- 


metry and the energy input only. Since it has 


now been shown [1, 8] that properly designed 
spray columns can be operated up to the flooding 
point with little or no coalescence of the dispersed 
phase, and that variations in droplet size over a 
wide range of nozzle velocities are not great [6), 
the holdup and flow rates are also related by 
equation (1). Consequently the droplet size and 
hence V, will be substantially constant for a 
given column and system, and therefore V, can 
be expressed as a function of the nozzle geometry 
and the physical properties of the system. 

If the flow rate of one of the phases is fixed, an 
increase in the flow rate of the other phase will 
result in an increased value of the holdup until 
It m iN 
be assumed that flooding occurs when V, or V. 


the flooding point is reached (see Fig. 1). 


reach a maximum value, L.e.. 
dV. da 0, 


curves beyond the flooding point are unrealizabl 


when dV, dr or 
Since those portions of the holdup 


in practice, the value of x corresponding to thi 
maximum flow rates also represents the limiting 
holdup at flooding, although this condition is not 
obtainable mathematically from equation (1). 


Carrying out the differentiations described 


above : 


7 dy) 0 


1) when (dV. 


Vea 2r,)(1 ry)? when dV drv)—0 (3) 


« 


The value of Py. the limiting holdup at the flooding 


point, is obtained by eliminating V, between 


equations (2) and (3) and solving for x, thus : 


(L? + sL)** — 31 


$(1 —L) 


“where L (4) 
A graphical solution of eq. (4) is given in Fig. 2. 

It will be seen from eq. (4) that the limiting 
holdup can be expressed in terms of the flow ratio 
only, and is independent of Vy. Furthermore, 
combining eqs. (2) and (3), and differentiating 
V, with respect to V.: 


dV. 

a aa 4 (5) 
dij c ] ay 
Consequently, the limiting holdup can also be 
evaluated directly from the slope of the flooding 


curve at any point. From the form of eq. (2) if 
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Graphical solution of equation (4). 


will be seen that a plot of V, versus #/ (1 ry) 
is a straight line through the origin of slope 2V’,. 
Such a plot is of use in correlating or extending 
incomplete flooding data, or conversely, if Vy is 
known for a particular system and nozzle con- 
figuration, it is possible to predict the entire 
flooding curve. Equations (2) and (3) can also 
be used to determine the limiting values of each 
flow rate as the other approaches zero. Thus 
when Vz 0, 2, > 0 and V, = V in the limit. 
- 0, V4 
Furthermore, it will be apparent 


Similarly, as V, a ~05 and Vy 
in the limit. 
from eq. (5) that when V, is plotted against V, 
as abscissa, the curve meets the abscissa tangen- 
tially as V, + 0, but meets the ordinate with a 
las V. +90, 


slope of These conclusions are 


borne out experimentally and will be discussed 
further in Section 4 below. 


(b) Calculation of V, 


The droplet sizes produced by different diameter 


nozzles have been studied fora number of aqueous- 
Hayworrn 
Treyeat [6], who obtained the equation : 


wd . 967 x 10 (tao) (") 
6 Ap 6 


194 >» 10 (72x) 
. Ap 


organic solvent systems by and 


+ 


= DE pow wen 2 
Ap Z 


which relates the droplet size to the nozzle 


+O 





diameter, physical properties of the system, and 
the velocity vy of the dispersed phase through 
the nozzle. 
Writing V, in terms of the droplet diameter [2] : 
PA 
V, = 5-36 x 10°98. (7) 
and TreyeaL evaluated k over a 
range of variables and obtained the ex- 


Hayworrn 
wide 
pression : 


k 


2-28 


10° Dg" 


vy? Oo 


oe (8) 


ff 
Substituting for & in eq. (7) and replacing d in 
terms of V, in eq. (6), it is seen that 


134 « 10° nad 


— Ap®? gi? vy? 


10 3 


,., 136 o-19 O75 4 
1 9.58 (“— Ke Dy —*) 


Ap*g 
y Dy 
Ap 
D a vy a7 ty -) 2 


Ap 


91 





10+ 

| 
Remembering from eqs. (2) and (3) that, at the 
flooding point when V, = 0, Vy = V,/4, then 


p — the 
a 


and by definition 


(10) 
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L, 


° 11 
aN (11) 


Combining eqs. (10) and (11), 


1 v 
taN 


substituting for (9), and 


2.29010 ‘| 4 | 
a. 
n D?™,* Pal A 
Ap*g 


10 (728) 0 
p 


} (12) 


x) | 


bt | 


‘» to be determined for 


haquation (12) enables J 


ny particular system, column and nozzle 
geometry, and is readily solved by setting the 
equation equal to y and plotting this quantity 
The value of V, 


solution of 


eainst assumed values of 1). 


corresponding to y 0 is then the 


3 APPLic Pusptisuep Dara 


ATION TO 
(1) is demonstrated 


% for the system ethylene 


The bility ol eq. 
Fig. 


dispersed into water, for which the 


applies 
graphically in 
dichloride 
oldup values were determined by the pressure 
drop method [8]. 

Some typical flooding curves for a number of 
aqueous-orga nic solvent s\ stems (1, 8] are shown 
in Fig. 4(a): these are extrapolated to zero flow 
The 


corresponding linear plots of eq. (2) for the same 


rates in accordance with eqs. (2) and (3). 


data are shown in Fig. 4(b), and in these cases, 
the values of vy employed were calculated from a 
knowledge of the flow ratios using eq. (4). All 
the available flooding data [1, 8] were plotted in 
the manner shown in Fig. 4(b), and the values 
of 4 
from the slopes of the resultant plots. A graphical 


, for the various systems were calculated 
comparison between these V, values and those 
caleulated from eq. (12) is shown in Fig. 5. 

\ graphical comparison between the limiting 
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Plot of experimental values of V, versus values 


calculated from eq. (12). 


holdup values determined experimentally and 


those calculated from eq. (4) is shown in Fig, 6. 


4. Discussion 


It will be 
determined both below and at the flooding point 


seen from Fig. 3 that holdup data 


can be related to the phase flow rates by eq. (1). 
There is some deviation from linearity below 8°, 
holdup, the reported values being on the low side, 
but this may well be associated with the difliculty 
of making eccurate measurements in this region. 
The slope of the plot gives a Vy value of 1250 ft, hr 
as opposed to a value of 1380 ft, hr calculated 
from eq. (12). If the flooding data only are plotted 
in accordance with eq. (2), as shown in Fig. 4 (b), 
the value of V, derived from the slope of the 
plot is 1200 ft/hr which is in good agreement 
with the value obtained from Fig. 3. 

The validity of eqs. (2), (3), (4) is well demon- 
strated in Figs. 4(a) and 4(b). In the first 
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Fic. 6. Plot of experimental values of ay versus values 


calculated from eq. (4). 


diagram some typical flooding curves are shown 
for a number of systems, and it will be seen that 
these extrapolate well to zero flow rates of each 
phase. Thus, as shown previously, when V, = 0, 
V,/4, and when V, = 0, V) = Vy, where 
the appropriate value of V, is obtained from the 
slope of the corresponding linear plot in Fig. 4 (b). 
A plot of the V, values obtained in this way 
against those calculated from eq. (12) is shown 
in Fig. 5. 

When determining values of V, from a plot 
of eq. (2), points corresponding to a dispersed 
phase velocity through the distributer nozzles 
greater than I1ft/see were discarded, since 
velocities in excess of this figure give rise to 
droplet sizes which are erratic and non-reprodu- 
cible [6]. It will be seen that the majority of the 
calculated values lie within + 20°, of the experi- 
mental values. Those values which lie outside 
this region tend to be on the high side, but this 
is to be expected, since Minarp and Jounson [8] 
defined flooding as the point at which the solvent 
phase was totally rejected from the column entry 
cone, rather than the true Elgin flooding point at 


V, = 


20 


which the solvent phase was just rejected from 
the column proper. 

Fig. 6 shows a comparison between values of 
the limiting holdup calculated from eq. (4) 
values obtained by the settling method, 


and 
and 
from pressure drop measurements across the 
column. There is good agreement between the 
calculated values of a, and those determined 


all 


cases the values determined by allowing the two 


from pressure drop messurements, but ‘n 


phases to settle out in the column are consistently 
These latter values, shown in Fig. 6, are 
by MrNArp 


high. 


values reported and JOHNSON, 


multiplied by the factor 0-92. This correction 
is necessary since, in the present case, we are 
considering the holdup in the column only, whilst 
the above workers defined their fractional holdup 
as the ratio of the height in the column occupied 
by the settled dispersed phase to the height of 
the interface above the dispersed phase nozzle [7], 
It is probable that holdup values determined by 
the pressure drop method are inherently more 
accurate than those determined by settling since, 
in the latter case, the dynamic level of the dis- 
persed phase in the disengaging section of the 
column was not amenable to accurate measure- 
ment. Consequently these values cannot be 
considered to have an accuracy greater than 

10%, [7]. 

In the cases of rotary annular and pulsed 
columns, the droplet size is a function only of 
the energy input, column geometry, and physical 
properties of the system, and is independent of 
the phase flow rates. In the present case, however, 
some dependence of droplet size upon the nozzle 
velocity of the dispersed phase has been demon- 
strated [6]. 
over a fairly wide range of nozzle velocities is not 


However, in general, the variation 


great, the droplet size tending to increase with 
velocity up to 10 cm ‘see and then decrease with 
Thus, over a tenfold 


velocity up to 30 em sec, 
range of nozzle velocities from 2 to 20 cm/sec, 
Haywortn and Treysar [6] reported droplet 
sizes for the water, benzene-alkaterge system 
varying from 9-3 to 10-1 mm using a 7-86 mm 
nozzle, and from 4-8 to 5-9 mm using a 1-55 mm 
nozzle. Since the nozzle velocity influences the 
droplet size by virtue of its contribution to the 
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kinetic energy of the issuing droplet, it seems 
reasonable to suppose that the influence of the 
continuous-phase velocity will also be small. 

With these considerations in mind, Vy was 
assumed to be substantially independent of phase 
flow rates in eq. (1), with the result that eqs. (2), 
(3), (4) are considerably simplified whilst repre- 
senting the data adequately, as shown in Figs. 4 
and 6. 

It is of interest to note that although there is 
a tendency for the droplet size to decrease with 
nozzle velocities higher than 10cm sec, the 
dispersed phase holdup increases, thereby in- 
creasing the probability of droplet coalescence. It 
may well be therefore that, at high solvent flow 
rates, these two opposing effects should result in 
a substantially stable droplet size distribution. 


5. CONCLUSIONS 


The dispersed phase holdup in a spray 
extraction column can be related to the super- 
ficial phase velocities and the characteristic 
droplet velocity by means of eq. (1). 


The holdup at flooding is independent of the 
droplet size, and can be predicted from a 
knowledge of the flow ratio alone by means 
of eq. (4). 


(c) Flooding data can be extended and correlated 
in terms of the characteristic velocity by 
means of eqs. (2) and (4). 


(d) 


The characteristic velocity for any given 
systern and nozzle geometry can be calculated 
with an accuracy of + 20%, from eq. (12). 
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NOTATION 


cross-sectional area of column - ft?. 
= cross-sectional area of dispersed phase nozzle - ft?. 
diameter of dispersed phase nozzle — ft. 
= acceleration due to gravity — ft /hr?. 
= number of dispersed phase nozzles. 
superficial phase velocity - ft /hr. 
characteristic velocity of droplets — ft/hr. 
= fractional holdup of dispersed phase. 
density — Ib /ft*. 
= density difference of phases — Ib /ft*. 
== viscosity — lb. /ft.hr. 
= interfacial tension of phases — Ib /hr?. 
fractional voidage of packing (unity for spray 
columns). 
Subscripts 
¢ = continuous phase 
d = dispersed phase 
f = at flooding 
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On the characteristics of a non-isothermal chemical reactort 
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Abstract—The temperature reaction history in a chemical reactor is analysed. A _ practical 
method is outlined for obtaining particular solutions of the pertinent differential equation for 
specific values of the parameters which control the heat transfer and the reaction kinetics of 
the system. With help of a phase plane analysis, properties of the general solution of the differ- 
ential equation are established which exhibit the influence of the above mentioned parameters 
on the temperature reaction history. The analysis indicates two physically distinct operational 
ranges of the reactor which depend essentially on the exterior heat transfer characteristics 
of the system. 


Résumé—L ‘auteur étudie mathématiquement lévolution de la température de réaction dans 
un réacteur tubulaire. I] ébauche une méthode pratique d’obtention des solutions particuliéres 
de l’équation différentielle appropriée pour des valeurs spécifiques des paramétres qui contrélent 
léchange thermique et la cinétique de la réaction du systéme. Par une analyse dans un plan, 
lauteur é¢tablit les propriétés de la solution générale de l'équation différentielle qui prouvent 
Vinfluence des paramétres mentionnés ci-dessus sur la variation de la température de la réaction. 
L’analyse montre l'existence dans le réacteur de deux domaines opérationnels physiquement 
distincts qui sont fonction essentiellement des caractéristiques d’échange thermique avec l’extérieur 


du systéme. 


INTRODUCTION 
Tue problem of determining the temperature 
history and the history of the reaction products 
in a flowing gas stream in which a chemical 
reaction occurs is of considerable importance in 
many fields of technology. It arises for example 
in the design of catalytic converters or homo- 
geneous gas reactors in the chemical process 
industry. The interest in the thermal aspects of 
continuous flow reactors stems from the fact that 
many reactions must be held within relatively 
narrow temperature limits to be feasible. This 
requires provisions for heat exchange to prevent 
the occurrence of temperature excesses or 
deficiencies. The history of the reaction products, 
on the other hand, indicates the degree of com- 
pletion to which the reaction has been carried 
out, and is thus an important parameter for an 
economic reactor design. The particular problem 
discussed below may be described as follows (see 
also [1]). An inert gas stream mixed with small 
quantities of reacting components capable of an 


exothermic reaction flows steadily through a long 
reactor tube of constant wall temperature. The 
entering gas stream which has been preheated to 
this temperature undergoes a chemical reaction 
as it passes through the tube. In its progress the 
mixed mean temperature of the gas rises at first 
with increasing distance from the entrance, since 
the amount of heat liberated by the chemical 
reaction exceeds the quantity of heat lost into 
the surrounding medium through the walls of the 
reactor. At a certain distance from the entrance 
the two effects attain an equilibrium and the gas 
temperature reaches a maximum. Thereafter, due 
to the consumption of reactive components, the 
mixed mean temperature of the gas decreases and 
attains the wall temperature on completion of 
the reaction. 

The primary purpose of this paper is to deter- 
mine the gas temperature as a function of the 
degree of completion of the reaction, the para- 
meters which control the heat transfer, and the 
chemical reaction characteristics of the system. 


: + This work was supported in part by the Office of Naval Research under Contract Nonr-222 (04). 
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Although the above description and the analysis 
given below is applied to exothermic reactions, it 
can with 
endothermic 
outlined is not new. A number of analyses have 
been proposed which account for both the 
chemical kinetics and the transport properties in 
a flowing gas stream [2, 3, 4, 5). 

In one model one assumes that the reacting 


minor modifications be extended to 


reactions as well. The problem 


gas possesses a uniform average velocity, con- 
centration and temperature distribution trans- 
verse to its direction of flow. In this one-dimen- 
sional theory the heat transfer is assumed to occur 
entirely across a fictitious thermal boundary layer 
at the wall of the reactor. 


transfer is described through the introduction of 


The process of heat 


an empirical transport coeflicient, the magnitude 
of which must be known a priori. 
detailed description, which is usually based on a 


In a more 


uniform velocity distribution, the temperature 
variation transverse to the direction of flow is 
taken into account. This avoids, then, the use 
of an overall heat transfer coeflicient. 

Either formulation of the problem leads to 
nonlinear differential equations of some com- 
plexity. Hence, nearly all solutions presented in 
the literature have been obtained by graphical 
or numerical methods for specific values of the 
parameters entering the problem. The results 
consequently lack the generality of an analytical 
solution. In this paper we shall be primarily 
interested in exhibiting both the qualitative and 
quantitative behaviour of the basic parameters, 
which characterize the heat transfer and reaction 
kinetics, on the average overall performance of 
the reactor. The model for the discussion will 
be the system described in terms of an empirical 
heat transfer coeflicient. This will, of course, 
leave the finer details dealing with the transport 
will be 


properties out of discussion. These 


presented in a subsequent paper. 


ANALYSIS 


The following assumptions are made. 

(1) Steady-state conditions exist in the reactor 
tube, which is of circular cross-section and con- 
stant wall temperature 7’, along its entire length. 

(2) The reactants form a slight admixture to 
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an otherwise inert gas stream so that the mole- 
cular weight of the gas passing through the 
reactor remains essentially unchanged. The gas 
as a whole obeys the perfect gas law. All gas 
properties with exception of the density and 
temperature are assumed to be constant through- 
out the reactor. 

(3) 
composition of the gas are assumed uniform at 
each cross-section of the tube, but vary in the 


The velocity, temperature, density, and 


direction of flow. The heat transport mechanism 
transverse to the direction of flow is taken into 
account in the form of Newton's law of cooling. 
Conduction and diffusion effects in the direction 
of motion are neglected. 

(4) 
gas and the reactor wall is relatively small, i.e., 
(T — T,,) < T,,. 

(5) The and 
assumed to follow the Arrhenius reaction rate 


The temperature difference between the 


reaction is monomolecular is 


law. Despite these restrictions the model des- 
cribed below exhibits the essential charactristies 
of a non-isothermal chemical reactor. 

In view of the above assumptions, conservation 


of mass requires that 
G (a constant) 


(1) 


where p is the density, u the velocity and G-the 


pu 


mass velocity of the gas. 
The kinetics of the reaction can be expressed 
in terms of « (0 < « < 1), the fractional comple- 


tion of the reaction, by 


de 


A(l 
dx ( 


Here 2 is the distance from the entrance of the 
reactor, A the frequency factor, E the energy of 
activation of the reaction, R the gas constant (in 
thermal units) and T the temperature of the gas. 
The expression for the conservation of energy 
for this model represents a balance between the 
heat carried by the gas stream, the heat liberated 
by the reaction, and the heat lost 
through the wall of the reactor tube, i.e. 


chemical 


de 


Qp, u de 


aa dT 
p dx 


4h 
~5(T-T.). (3) 


Here Q represents the heat of reaction per unit 
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mass of the reactant converted. p, is the initial 
density of the reactant, A the overall heat transfer 
D the diameter of the reactor tube 
and ¢ the specific heat of the gas. 


coeflicient, 


In addition one requires the equation of state 


PM | 


4 
R, |} T *) 


p 
where M is the average molecular weight, P the 
gas pressure, and R, the universal gas constant. 
The associated initial conditions for equations (2) 
and (3) are 


«(0)=0; T(0) = T,,. (5) 
The equations (1) to (5) define the model with 
which we shall be concerned. It is a mathematical 
model from which one can deduce the influence of 
the most important physical factors on the 
progress of the chemical reaction. 


In view of assumption (4) the temperature rise 


in the gas stream is moderate enough so that one 


if, £-f), 
| A | ( i) 


can write 


The dimensonless excess temperature 
0 = (7) 


having been introduced, the equations (2) to (4) 
become 


Ae ®/¥Tw (1 ee? (8) 


th RTZ 
DE 


PM 
RoT 


0 (9) 


RT. @) where p, 


10) 
E ( 


The independent variable x can be eliminated 
between equations (8) and (9) with the result 

: RT2Z d@ Q th RTZ 1 

Ede “ D E Ac®®Tw\ 


de 6 
ve > (11) 
€ 


On eliminating the density p with the help of 
(10) one obtains finally 


dé 


where 


¢ , 
| . Pr) ( : Jand 8 th I : . (13) 
CT pe) \ RT. De A pye © BT « 


Equation (12), which is best suited for our 
purposes, will when integrated yield the dimen- 
sionless gas temperature @ as a function of the 
reaction history ¢« in terms of the essential 
parameters x, 8 and RT, E which describe the 
heat transfer and reaction characteristics of the 
system. The initial condition associated with 
equation (12) is, in view of (5), 


0 (0) 0 (14) 


A general analytical solution of equation (12) in 
terms of elementary functions can unfortunately 
not be given. It is possible to obtain a series 
expansion in the parameters x, 8 and RT, E but 
little can be learned in this way about the essential 
nature of the solution. A more useful approach 
proceeds by a geometric phase plane discussion 
which, incidentally, yields also a_ practical 
method for solving equation (12) for any specific 
values of the parameters. 

For the construction of a particular solution 
curve of equations (12) and (14) one can employ 
the isocline method. A lineal element field of 
equation (12) is constructed from 

I, B _ se K (aconstant). (15) 
(1 RT, 7 "1l—ef 
E 


Solving this for « yields a family of curves, the 
isoclines 
B/« 
; (x _ RT, 0 
E 


€ be* 


(16) 


which are plotted in a @, ¢« plane for a specific 
value of 8 2, RT, E, x, and for different values 
of the parameter K. An isocline for a specific 
value of the parameter K is the loci of lineal 
elements of constant slope K, and the solution 
curve must cross this isocline with this slope. 
Such a field is shown in Fig. 1 for the following 
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Fig. 1. 
equation (12) for 8 « 


The isocline field and particu'ar solution curve of 


3,2 = 10, RT,,/E = 1/10. 


values of the parameters: 8/2 = 3, « 10, 
RT, E=1,10. The solution curve, starting 
from the initial point (0-0) with a slope of 10 
according to equation (12), rapidly attains a 
maximum, and falls subsequently to end in the 


point (0-1) as required by the completion of the 


reaction. By this graphical solution method any 
particular solution curve can be readily con- 
structed and the temperature reaction history be 
determined. The construction of the isoclinal 
field becomes particularly simple if the term 
RT,0 E can be respect to 
unity, since in this case the members of the 
family (16) differ merely by a scale factor. 
Having exhibited this convenient method for 


neglected with 


determining particular solutions for specific values 
of the parameters, we turn next to a detailed 
discussion of the properties of the general solution. 
Consider for this reason the complete isoclinal 
field of the differential equation (12) as shown in 
Fig. 2 for some particular values of the para- 
meters a, 8 a and RT, E. Within the circular 
sectors the isoclines are the loci of lineal elements 
of either positive or negative slope K, as indicated. 
Starting from the top and to the left of the line 
« = 1 and moving in a counterclockwise direction, 
the value of K increases from minus infinity to 
zero, changes sign and then increases from zero 
to plus infinity at the bottom of the line « = 1. A 
similar behaviour is encountered to the right of 


« = 1. All isoclines intersect at the point (0-1), 


ry) 0 GOCLINE 








Fic. 2. The phase plane of equation (12). 


which is a singular point of the differential 
equation (12). 

The first question which comes to mind con- 
cerns the structure of any particular solution 
curve originating from the ordinary point (0-0) 
with slope « in accordance with equations (12) 
and (14). Will the solution curve terminate in 
the required point (0-1) or is it possible for it to 
miss this point? In the latter case only three 
distinct possibilities can arise as shown by curves 
A, B, and C in Fig. 3. In case A assume that the 
solution curve meets the ¢ axis somewhere between 
In order to do this it must cross 
would have a 


zero and one. 
the K 
minimum or an inflection point with zero tangent 
at the crossing point. In the first case the curve 


= 0 isocline, and hence 


would turn up, which is clearly impossible since 
the slope of the lineal elements above the isocline 
is negative. In the second case the solution curve 
would continue with a negative slope into a 
region where the lineal elements are positive. This 
rules out case A. 

In case B, assume that the solution curve does 
not intersect the required point (0-1) because the 
curve turns upward and becomes asymptotic to 
the line « = &, (0 <€ <1). This is clearly 
impossible, since the slope of the solution curve 
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INCREASING K 
KO SOCLINE 


6 











€=! 
Fic. 8. Properties of the solution field in the neighbour- 
hood of the singular point. 


would be positive in a region where it must be 
negative. The third possibility is that the solution 
curve C assertedly intersects the line « = 1 at a 
point above the singular point. However, this 
is impossible, since the solution curve on approach- 
ing the line « = 1 from the left takes on an infinite 
negative slope, and hence cannot cross in any 
continuous fashion to the other side where the 
slope is positive and infinite. Consequently, the 
solution curve on approaching « = 1 will turn 
downward with infinite slope. Furthermore it 
cannot turn away, as shown by the segment D. 
Hence it must lead into the singular point. How- 
ever, the solution curve with infinite negative 
tangent is not the only possibility of reaching 
(0-1) as we shall see by investigating the specific 
nature of the singularity. One concludes from the 
above that, for all finite values of the parameters, 
equation (12) has a solution which, starting from 
the point (0-0), must end in the singular point (0-1) 
To gain further insight into the properties of 
the solution curves in the neighbourhood of the 
singular point, consider the appropriate form 
equation (12) takes close to this point. One can 
show in general [6] that the differential equation 


dy _ax+ by + P(2, y) 


= z (17) 
dz ca + dy + Q(2,y) 


in which the constants a, b, c, d are such that 
ad — be + 0, and in which P and Q vanish as 
x +y* does when z,y +0, has as its only 
singularities at the origin (2 = y = 0) those of the 
much simpler equation 


(18) 


Equation (12) can be brought into this form by 
expanding the exponential term about 6 = 0 and 
neglecting all terms higher than the constant 
term. In addition, the co-ordinate origin is 
shifted to (0-1) by introducing the new variable 
(1 —e). Hence 


d0 BO —a(1 —«), (19) 
(1 — e) 


d(l—e) 
Without going into the details of the analysis 
here, it can be shown that the point (0-0) of 
equation (19) is a nodal point which has the 
characteristic that all solution curves passing 
through the immediate neighbourhood of this 
point will end in the singular point, as already 
found above. The shape of the solution curves 
near the singular point, aside from the trivial 
characteristic « = 1, is readily established by 
integrating equation (19), with the result 
, B x 
Ae) = C(l —€) 234" e), B +1 (20) 


We) = C(1 —e) —a(1 — 6) In(1 — 6), B=1. 


From equation (20) one observes that for 8 > 1 
the solution curve @(«) becomes tangent to the 


)as e > 1. On the 


-'@ 
B-—1 

other hand if 8 <1, the equations reveal that 
the slope of the solution curve ending in the 
singular point becomes negative and infinite. 
These properties allow one to define two distinct 
operational ranges of the reactor in terms of the 


straight lines of slope 


parameter £. 
Equation (13) shows that 
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is the quotient of two reciprocal times. 7, is 


a characteristic relaxation time for the cooling 
of the reactor bed due to the overall heat transfer 
mechanism through the walls of the reactor, 


while t.4.~ iS a measure of the speed of the 


chemical reaction. A small value of 8 implies a 


small value for and a large value of + 


‘chem 


or. in other words, a fast chemical reaction and 


conv 


a slow heat transfer exchange with the outside. As 
shown above, for values of 8 < 1, the solution 
curve ends in the singular point with infinite 
negative slope. This implies that the temperature 
of the gas falls primarily as a result of the final 
> 1, 


and not so much as a result of the overall heat 


consumption of the reacting material as 


It seems appropriate to designate the 
the tem- 
perature maximum occurs close to ¢ 1, 
the 
In this regime the temperature 


loss. 
) 


operational range 0 < 8 <1, where 
as the 
range in which reactor behaves almost 
adiabatically. 
history is primarily determined by the progress 
of the chemical reaction and to a lesser extent 
by the overall heat transfer characteristics of the 
system. 

On the other hand the solution curves for the 


» 


operational range 1 8 < @ enter the singular 


point with slope (, =}, which can result in 


solution curves of varying degrees of steepness 
with the 
anywhere in the range 0 < « < 1. 


temperature maximum occurring 
This indicates 
that for this regime the exterior heat transfer 
plays an important role and is apparently the 
norma! for most industrial 
reactors. 

The two cases are exhibited in a phase plane 


Here the A 


shown 


operating range 


0 isoclines of 
for 


diagram in Fig. 4. 
(12) 


equation are schematically 


, AR phy 
DE 


> ~~ —-— —---- 
Ae ® *¥e 5 ~ Heat generation per unit volume by the chemical reaction 


while 


Q p, 
© Pw 
RT? 
E 


Heat loss per unit volume to the outside 
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€=! 


Fic. 4. The K = 0 isoclines for different values of 8 /«. 


different ratios of 8 x. 
the operating condition 8B = 4,2 = 1, RTy E 

1/10. The solution curve labelled A leaves the 
origin with slope 1, crosses the isocline Ba = } 
1 with zero slope, and descends 
steeply to the singular point as shown. 
trast to this almost adiabatic behaviour, consider 
the reactor in the normal operating range 8 = 120, 
a = 15, RT,/E = 1/10 shown by curve B. The 
increased heat transfer decreases the temperature 
maximum considerably, and shifts the latter to 


Consider, for example, 


close to «€ 


In con- 


a much smaller value of « than in case A. 
From Fig. 4 one can furthermore deduce the 
qualitative dependence of the solution on the 


Equation (13) shows that 


parameters a and ~- 
a 


these parameters are a measure of 


(22) 


Adiabatic temperature rise in the gas 
Effective temperature rise in the gas ~ 
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Consider a fixed « and two different values of 
(B/a). The two solution curves labelled C and D 
are shown in Fig. 4. In agreement with physical 
expectation the solution with a smaller value of 
(8/«) attains a larger temperature maximum. One 
verifies, on the other hand, for a fixed value of 
(B/a) and two different values of «, that the 
solution curve with the larger « will have the 
larger temperature maximum. If one considers 
the maximum value of @ as a function of « and 


B/a, i.e., @ 


2A max 


and 
52 


a a 
parameter which is held constant. The main 
conclusion in this case is that the temperature 
maximum in the reactor can be decreased with 
decreasing « and increasing 8/«. One can sub- 
stantiate this in the following limiting case which 
yields the temperature maximum as a function 
of the parameters. It is seen from Fig. 4 that for 
the complete range of « the maximum temperature 
is essentially controlled by the K = 0 isocline 
which depends only on the value of 8/x. If the 
isocline for K = 0 intersects the line « = 0 at, 
say 0, then clearly 0... < 9% for all « One can 
improve this estimate by the following con- 


= 0 (a, B/x), one has (mex) y <0 


, Of Sa 
<0, the subscript indicating the 


siderations. 

For a sufficiently large value of 8 « the maxi- 
mum reactor temperature will be small enough 
for one to neglect the term RTy@/E with 
respect to unity in equation (12). Physically this 
amounts to the assumption that the density and 
velocity variations through the 
negligible. Equation (12) can then be written 


reactor are 


1 dé 
(;)“ 


However, for this asymptotic condition one may 
neglect the derivative term on the left hand side 
almost everywhere except in the neighbourhood 
of the origin. Hence, for the greater part of the 
reaction history the 0, « relation is given by 


B Bg 


(24) 


e=} b6e* ~1— 


since @ is everywhere small. This is essentially 
the equation of the K = 0 isocline, which the 
solution curve follows closely. In the neighbour- 
hood of the origin the solution of equation (24) 
will have the form 


0 = ae + 0(e*). (26) 


The value of @ at the intersection of the two 


curves approximates @,,, and is given by 


This asymptotic relation verifies the dependence 


of the two parameters on @ as 


= discussed 


above. 


NOTATION 
frequency factor of monomolecular reaction, 1/h. 
specific heat of gas, BTU /Ib. °R. 
reactor tube diameter, ft. 
energy of activation, BTU/mol. 
mass velocity of gas, Ib. /sq. ft. h. 
heat transfer coefficient, BTU /sq. ft. h. °R. 


slope values of lineal elements, defined by equation 
(15), dimensionless. 

average molecular weight of gas, lb. /mol. 

gas pressure, lb. /sq. ft. 

heat of reaction per unit mass of reactant converted, 
BTU Ib. 

gas constant, BTU /mol. °R. 

universal gas constant, ft. lb./mol. °R. 

gas temperature, °R. 

reactor wall temperature, “R. 

time, h. 

gas velocity, ft./h. 

distance from entrance of reactor, ft. 

parameters defined by equations (13) or (22) and 
(23), dimensionless. 

degree of completion of the reaction, dimensionless | 
excess gas temperature defined by equation (7), 
Omax its maximum value in the reactor, dimension- 
less. 

density of gas, p,, its value based on reactor tube 
temperature, Ib. /cu. ft. 


, = initial density of reactant, Ib. /cu. ft. 





Paut L, Cuamprt 


REFERENCES 
Cuamere, P. L. and Grossman, L. M.; Appl. Sci. Res. Sec. A, 1955 5 245. 
Wiutson, K. B.; Trans. Inst. Chem. Engrs. (London) 1946 24 77. 
Grossman, L. M.; Trans. Amer. Inst. Chem. Engrs. 1946 42 535. 


Sincer, E. and Wirneiom, R. H.; Chem. Eng. Progress 1950 46 343. 


Baron, Tu.; Chem, Eng. Progress 1952 48 118. 


Sroxer, J. J.; Nonlinear Vibrations, Interscience, New York, 1950, Ch. IIT. 





Chemical Engineering Science, 1956, Vol. 5, pp. 217 to 225. Pergamon Press Ltd, 


On the condensation of a vapour of a binary mixture in a vertical tube 
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Abstract— General differential equations for the condensation of a binary mixture in a cylindrical 
condenser are given. The solution for a specific case, viz in which the main body of the vapour 
is in turbulent flow, is presented. For simplified experimental conditions an equation presenting 
the temperature difference between the main vapour stream and the wall, as a function e.g. 
of the amount of condensate and of the partial coefficient of mass transfer in the vapour phase, 
has been derived. 


Experiments on vapour mixtures of benzene and toluene, in a vertical condenser cooled 
with a liquid, have confirmed that the theoretical approach given can in some cases be used 
for the calculation of the condensation of a binary mixture. It also appeared that the partial 
coefficient of mass transfer in the vapour phase, which could be calculated with the aid of the 
equation derived, above a Reynolds number of 2,000, was a linear function of Re’. 


Résumé—On donne les équations différentielles générales pour la condensation d'un mélange 
binaire dans un condenseur cylindrique. On trouve la solution pour un cas spécial, a savoir 
le cas of: la masse principale de vapeur est mise en courant tourbillonnant. Pour des circonstances 
expérimentales simplifies, on a dérivé une équation dans laquelle la différence de température 
entre le courant principal de vapeur et la paroi est représentée comme une fonction p.e. de la 
quantité du condensat et du coefficient partiel du transfert de la masse dans la phase de vapeur. 

Des expériences effectuées sur des mélanges de vapeur benzéne-toluéne dans un condenseur 
vertical, refroidi par liquide, ont confirmé que l'approximation théorique donnée peut, dans 
quelques cas, étre employée pour calculer la condensation d'un mélange binaire. On a trouvé 
également que le coefficient partiel du transfert de la masse dans la phase de vapeur qui a pu 
étre calculé a Taide de léquation dérivée, était, au-desus du nombre Reynolds de 2,000, une 
fonction linéaire de Re®®, 


Tue knowledge on the condensation of vapours’ For the condensation in a vertical tube of a 
of binary mixtures is of importance not only in diameter, large compared with the thickness of 
itself but also for every non-adiabatic distillation the liquid film, the following differential equation 
of such mixtures. Only a few articles on this can be derived (see Fig. 1): 


subject have appeared [1, 2, 3, 4]. From these ° 
the approach of CoLpurn and Drew [4] to the 
present authors seemed interesting and will be 


included in the considerations put forward in 
the following. 


THEORETICAL CONSIDERATIONS 


One Component 





For the film type condensation of a vapour 
consisting of a single component with the liquid 
film in laminar flow, the well-known Nussevr [5] 
equation, existing in several forms, can be used. Fic, 1. Condensation of an unair system. 
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dL 
(AN)r, S dz 
After substituting: 6, . in which 
r = ‘ 
ps 
one obtains : 
By M\'® (AH , 
(Tp Ty) dz | * “4 (0 = rp L' . aL. 
p~ gs Ss A 


mean 


(2) 


difference 


AT 


temperature 
the wall: 


Introducing a 
between the vapour and 
(Tp 7 T wm defined by : 
Zz 
l . 
(T 
Z. D 


0 


AT Ty ) dz 


and taking for », p and A constant mean values 
at the average film temperature, integration of 
(2) between z = Z and L = Ly and z = 0 and 
L = 0 yields : 


1 ’ ead) J L; ’ 


S* 3 p*g 


1 (All)7), 
1-082{ ~ }- 
AT O8 (5) ; 


or: 


AT = A.L,“* 


Two Components 


Consider a vapour of a binary mixture which 
is partly condensed by flowing this mixture 
through a cooled vertical tube. The following 
assumptions will be made : (1) film condensation 
occurs ; (2) the liquid is completely miscible ; (3) 
the liquid is in laminar flow; (4) at the interface 
liquid-vapour the deviation from equilibrium 
is neglible ; (5) the total molal concentration in 
each phase is constant ; (6) the molal specific heat 
C, is constant; )7) stationary conditions occur, 

Two cases, viz. the vapour in laminar and 
turbulent flow, will be considered. 


Vapour in laminar flow 


Considering the vapour-phase and assuming 
that the molal flux (number of moles passing 
through the unit of area per unit of time) in the 
vertical direction (co-ordinate z) is p and in 
the radial direction (co-ordinate r) is n, the follow- 
ing differential equations can be derived : 


HreERTJES 
From the total material balance 


: (5) 
~ 4 
From the material balance for the more volatile 
component A, 


A solution for these equations has not yet been 


found. It may be noted that the equations 


(6) and (7) are identical if the ratio 1; 


A 
C,DC 


that is, if the Lewis number 1, because 


A 
Cp pD 
7 + p- 


Vapour in turbulent flow 


If the greatest part of the vapour is in turbulent 
flow, and if it is assumed that in this region no 
horizontal temperature or concentration gradients 
exist, the equations derived can be applied to 
the small layer flowing in viscous flow around 
the turbulent kernel. (See Fig. 2). 

















¢ 

ow 

Fic. 2. Condensation of a binair system with the main 
body of the vapour in turbulent flow. 
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The thickness of this layer related to the 
diffusion is 5p, related to the heat transport 5,y. 
It may be remarked that 5,/5, = (Le)'”*. 
Thus for Le = 1, dy = 8p. 

When this thickness is small compared with R, 
and assuming that p is constant in vertical 

on 


direction, equation (5) yields : — = @, 
or 


and 


thus n is constant (= N) in radial direction. 
Moreover, for the thin layer, the temperature and 
concentration gradients in vertical direction 
can be neglected, compared with those in radial 
direction. Therefore the equations (6) and (7) 
can be simplified to 


(8) 


0. (9) 


When the diffusion in the liquid film is negligible, 
the composition of the formed condensate is 
the same as the composition of the liquid at the 
the interface, 7g. Observing that (if r is taken 
from the “ boundary ” of the turbulent kernel) 
for r= 0, y 5p. ¥ = Yq. We 
can derive from equation (8) : 


N 
DC 5p 
€ . 


Yp. and for r 


% — YG 
%— Up 


(10) 


From equation (9) it follows, because for r = 0’ 


T = Tp, and for r = dy, T = Tg, that 


(11) 


and by condensation (q,) is 
qe = N (AH)p,. (13) 


Therefore the total heat flux through the inter- 
face vapour-liquid (q,,,) is 


Got = Va T+ Ye (14) 


Equation (10), (12), and (13) have already been 
given by Co_purn and Drew. Expanding the 


D 


above by rearranging and introducing k = 5.” 
D 


equation (10) yields 


(15) 


If partial condensation is considered over a not too 
great range, the vapour line in the 7-2 diagram 
may be assumed to be straight, and if 7” is the 
temperature of saturated vapour of a composi- 
tion equal to ag (see Fig. 3). 
Yo — Ya _ Tp — T; P 


ae (16) 
t—Yyq TT’ — Tg 


From equations (15) and (16) follows after 
rearranging, 


N 


Typ —Tg=(T' - T) (1 - e*). (17) 


Substituting equation (17) in equation (12) 


A 
and introducing 5. h, one obtains, because 


w 
the Lewis number is taken, unity, and therefore 


-(T’ —T,) NC, (18) 


and therefore 


N(AM)z, + NC,(T’ — Te). (19) 


Trot 
(AH); is the enthalpy difference between 
saturated liquid and vapour at equilibrium at 
the temperature 7g. Introducing an enthalpy 
difference (AH Jag (see Fig. 3) between saturated 
liquid (7g, ag) and the saturated vapour of the 
same composition (7", zg), equation (19) sim- 


plifies to 
Gur = N (AHN), (20) 
because (see Fig. 3): (AH)7, 


+ C,(T" — T<) —_ (AH), 
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‘1G. 3. Schematic representation of the enthalpy and 
boiling point as a function of the composition. 


and it has been assumed, that C,, is constant over 
the temperature range involved. 

For the condensation of a binary mixture the 
height dz, 


therefore 


differential equation over a_ short 


under the conditions assumed, 


be written (see eq. 2) 
| 3y M ) * (AH) x¢ 


p2-S-g S-A 


may 


L’? dL, (21) 
Introducing again a mean temperature difference 
between interface and wall defined by 
Zz 
1 ” 
Ta.— 7 ln = | (7. Ty) dz, 
(46 " Z. G " 
v0 
and taking mean and constant values over the 
whole liquid film for p, », A and M and over 
the surface for (AH),_., integration of (21) yields 


(Te — Tw)m = 4,47" (22) 
» M 


1 1/3 
1,082 (Al/);,.- Zz (". sys] : 


in which A, = 


For the vapour it has been found that : 
N 


Tp — Tg =(T’ — To) ( —e*). (17) 


M. Heertses 


Because only partial condensation will be 
considered over a rather small range of con- 
centrations, (7 — Tg) can be taken as fairly 
constant over the whole pipe. The same holds 
for k. The amount of vapour condensing N will 
be determined by q,,,. which in turn will be 
determined by the overall coefficient of heat 
transfer for the whole apparatus and the tem- 
derature difference between vapour and cooling 
fluid, 
described below, a liquid was used as a cooling 
fluid. 
liquid with that in the condensing vapour, it 
can be concluded that bere the greatest resistance 


For the condensation experiments to be 


Comparing the resistance in the cooling 


against heat transfer was in the cooling liquid. 
Therefore the heat transfer coeflicient 
will but slowly decrease from the top of the pipe 
to the Moreover, the temperature 
difference between liquid and vapour was increas- 


overall 
underside. 


ing somewhat from top to bottom, therefore the 
product of overall coeflicient and temperature 
difference, that is q,,,. can in approximation be 
Therefore N is constant 
Lz 
S.Z 
Tg is 


taken as constant. 


over the whole tube and equal to N 
This leads to the conclusion that Tp 
constant for the whole tube. 

The 
vapour and wall AT 


difference between 


T yw)» therefore is 


mean temperature 

al 

(Tp 
Sn 


AT = (T" oo) + A, Lf/* (28) 


T;;) (1 


Assuming that the mean composition of the 
liquid film at height z is 2, a material balance 
over the element dz gives 


d (aL) = «gal. 


As L = NSz, because under the condition 
assumed N is constant over the whole length 
of the tube, the equation gives 
d (az) = &@ dz, 
Zz 
from which follows az . Z | % dz. 
0 
Zz 
. = 1 ; 
As %@ can be defined by %% Z fr dz, it follows 
that “¢ 
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= Vz. (24) 


In other words, the mean composition of the 
liquid film at its surface is equal to the mean 
composition of the liquid leaving the apparatus, 
If it is wanted to check the relation found in 
equation (23), this can be most easily done at 
constant %, (or what is the same, at constant 2) 
and therefore constant (7" — Tg.) 

As the ratio -5 Es 


e eq. 22) for benzene 
7—= 3 (see eq. 22) for benzene 
p.S”".g. A 
and for toluene is practically the same, it has 


assumed that this is also true for their 
mixtures. Therefore A, is only a function of 


(All), Under the the 


enthalpy as a function of the composition is a 


been 


z 


assumptions made, 


straight line. Therefore 


A, = A4tg + Ay(l — 4%). 


At constant @,, also A, is therefore constant. 
The curve of AT versus L} ® (eq. 23) at constant 
@, (or ®z) will therefore depend only on k, because 
for a fixed pressure C is constant. For k = « 
the system will behave as an unair system. 
Equation (23) then simplifies to 

AT = A, . L#", 


" (26) 
the Nusselt equation, as could be expected. 
For Lz © — 0, equation (25) simplifies 


to 


oO,ase 


AT = T’ —T, + A, Lz". (27) 


All the curves for values of k between 0 and 
will have equation (27) as an asymptote. In 
Fig. 4 the form of several curves for a fixed 2% 


and therefore 7’ — Tg has been drawn. 


at 








Lz 1 


Ha(T-Tgt-e 52 Key say Lgl 


Fic. 4. Schematic representation of the equations (23), 
(26) and (27) at constant zg. 


EXPERIMENTAL 
To check the above considerations, the con- 
densation of a mixture of benzene and toluene 
carried out. The principle of the 


has been 


measurement was as follows (see Fig. 5). 
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The apparatus, 


A vapour was generated in a cylindrical boiler by 
electrical heating of the liquid and introduced at the 
bottom of the inner cylindrical tube (diameter: 2 R: 
88 mm, height Z : 1080 mm) of a double pipe condenser. 
The wall of the inner tube was 3mm thick. The vapour 
was partially condensed by means of a cooling liquid 
(a lubricating oil of low viscosity) introduced tangentially 
in the outer pipe. The vapour which did not condense 
on the inner pipe was liquified in a total condenser cooled 
with water. The liquid formed here was collected, its 
quantity measured, and returned to the boiler. The 
condensate formed in the experimental condenser was 
collected, brought outside the apparatus, its quantity 
measured, and returned to the boiler. The cooling liquid 
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was pumped through the apparatus and was kept at a 
constant temperature by passing it through a thermostat. 
Eight copper-constantane thermocouples were installed 
in the apparatus, four in the centre of the inner condenser 
tube (vapourcouples), four in the wall of this condenser 
tube (wall couples). The e.m.f. of these thermocouples 
was measured with a Dieselhorst compensator with five 
decades. The accuracy of the measurement was 0, 5 » V 
(corresponding with about 0-02°C). 

Samples of the vapour entering the experimental 
condenser could be drawn via a side tube by total con- 
densation. The pipeline connecting the two condensers 
was isolated, but no electrical heating to prevent con- 
densation was applied. The formed here 
was collected before it could enter the experimental 
condenser and was added to the product of the total 
condenser. 

At intervals samples were drawn from the streams 
entering and leaving the experimental condenser, and 
After analysis the samples were 


condensate 


these were analysed. 
returned to the boiler. 

A careful calibration as well as recalibration at intervals 
of a few months of the thermocouples used was necessary. 
This calibration has been carried out as follows. 

The vapour couples, which were easily detachable from 
the apparatus, were calibrated in the vapours of pure 
boiling liquids. The liquids and the constants used were: 





boiling point °C | dT/dP, “€ 


at 76cm Hg 


mm Hyg. 





56-131 0 0386 

64-65 0-035 

80-122 00-0426 
100-000 0-037 
110-606 O-O46 


acetone 
methanol 
benzene 
water 
toluene 











The fixed wall-couples were calibrated in the apparatus 
on the vapour couples. For this purpose and 
total condenser were detached from the experimental 
condenser, the ends of which were closed with a flange. 
The whole condenser was very heavily isolated, whilst 
the top was also heated electrically to compensate for 
heat loss through the beams from which the apparatus 
hung. During at least 5h, oil of a constant temperature 
was pumped through the annulus of the experimental 
condenser until the vapour couples all had constant 
and the same temperature (deviation of 0-05°C). The 
mean of these was taken as the actual temperature. 
The wall couples were read at the same time, and reading 
was continued over a rather long period. At least 20 
readings were taken for each couple. The mean of these 
with a maximum deviation of 3 nV was taken, and assumed 
to be the same temperature as the vapour couples. The 
calibration of the wall couples was carried out at the 


boiler 
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temperatures 70, 80, 95 and 110°C, The calibrations 
were in every case evaluated in lines of the form 


E =aT + bT?. 


In the apparatus the condensation of benzene and 
toluene and their mixtures has been measured. The 
compositions of the binary mixtures were measured by 
means of a semi-micro boiling point apparatus. The 
accuracy was 0-2%,. 


and toluene are 


AT in °C and L, 


The results for benzene 
represented in Fig. 6 with 
in moles /sec. 
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Fic. 6, The results of the condensation experiments with 
benzene and toluene. 


Excluding the observations at very low amounts 
of condensate the measurements can be correlated 
by the equations 


AT = 0-4 
AT = 0-5 4 


220 L,** for benzene, 
340 L,*” for toluene. 


The straight lines do not pass through the 
origin as is demanded by theory. The cause of 
this seems to be that the amount of thermocouples 
used is too small, and moreover that the wall 
couples are not all well placed, causing a changing 
and too high a temperature difference for the 
larger amounts of condensate. 

The theoretical slope of the lines can be 
calculated with the aid of equation (3). 

The factor A has been calculated from the 
physical constants available and determined for 
benzene at 79°C, for toluene at 109°C, because 
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the temperature difference between wall and 
vapour is 8°C at a maximum. 
The following table gives the constants used. 





Benzene Toluene 





7370 7960 
3-5 10+ 3-34 
325 10% 2-55 
0-815 0-783 
78 92 
185 210 


cal /mole 
10+ cal/°C cm.sec 
10°% g/cm. sec 
g/cm® 

gz/mole 

°C /(moles /sec)* # 








As can be seen, the differences in slope between 
the calculated and determined lines are 
appreciable. This may be due to inaccuracies 
in the constants used and to deviations from 
physical behaviour assumed. 
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(Lz in molesisec) 
Fic. 7. The results of the condensation experiments with 
a mixture of benzene and toluene (series I, 7z = 0-33). 
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The experimental values as found will be 
utilized to evaluate the results of the condensation 
of the binary mixtures for which the same benzene 
and toluene has been used, and therefore will be 
used as a calibration of the apparatus. For the 
binary mixture two sets of experiments has 
been carried out, the first with a vapour contain- 
ing 0-47 mol fraction of benzene, the second with 
a vapour containing 0-75 mol fraction of benzene. 
The vapour velocities and the temperature of 
the cooling liquid were varied. 

In the first series, % ranged from 0-490 to 
0-268 (mean value 0-33). The mean vapour 
velocity G (in moles /sec) ranged from G. 10* = 4 
to 87, the Reynolds number of the vapour from 
530 to 12070. In the second series %z ranged 
from 0-692 — 0-0,551 (mean value 0-61), G. 10° 
from 18 — 92 and the Reynolds number of the 
vapour from 1050 to 12360. 
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¢ 4 
———————> L75x10 
(Lz in moles/sec) 


Idem (series II, zz = 0-61). 





d. P. 
Although the values for 2, differ rather con- 
siderably, also in each series, the differences in 
A, (280 to 307 : 300 for series I and 256 
to 274: mean 260, for series IT) as well as in 
T’ — Tg, are of such a magnitude that a com- 
parison of the values obtained 
series is permitted. The 
plotted in Fig. 7 for series I, in Fig. 8 for series IL. 
The values for about equal Reynolds numbers 


mean 


within each 


results have been 


marked in the same and have 
The 
caleulated for @, = 0-33 and 
and the corresponding 7” 


the vapour-liquid equilibria data for benzene- 


have been way 


drawn have been 


0-61 re spective ly 


been combined. lines 


- Tc as read from 
toluene. The temperature differences are given 
in the figures as they have been found experi- 
mentally. The equations given in the figures for 
k © and k 


deviation in AT of the pure substances (see the 


0 are corrected for the observed 


equations above). 

As is clear from Figs. 7 and 8, the qualitative 
picture obtained is in agreement with Fig. 4 
given in the theoretical part. 

From the data collected it is possible to calculate 
Tw), Tz. 
and Ly, have been measured ; moreover A, can 
be calculated from equation (27) and 7” — T, 
read from the equilibrium curve. It can be 
found that : 


k for each experiment, because (7'p 


L,/S .Z2.C 
— "ny s/n (28) 
T,, — Al Ls 





' @ =seriesl 
x =seriesD 


viscous ~+tturbulent + 
flowm— - t —fiew-— 











Fic. 9. The partial coefficient of mass transfer (4) as a 
function of Re’ 5, 
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The form under the In sign in equation (28) is 
the ratio of the vertical distance of the point 
plotted to the asymptote, and the vertical distance 
between this asymptote and the line for k Dn. 
Therefore & is independent of the corrections 
mentioned before. 

The values found for k have been plotted in 
lig. 9 versus the Reynolds number of the vapour 
stream to the power 0-8. For turbulent flow 
in a pipe Sh: (Re)’* (Sc)'/*, If it 
that D and Se'* are constant for the systems 
investigated, a straight line will be obtained in 


is assumed 


such a case. 
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series 
series. 
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0.75 0,25 ° 


x7 ——_____—_ 


0,50 


Fic. 10, The temperature 7, calculated from T= T yw + 
A, L;** as a function of zz. 
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As can be seen from Fig. 9, this relation for 
the region of Re > 2000 (Re°* > 450) is fairly 
well satisfied. 

Finally 7, has been calculated from the data 
obtained and plotted versus #7 in Fig. 10, together 
with the lines for saturated liquids and vapours, 
Although most of the points lie below the line 
for the saturated liquids, the correlation is not 
too unsatisfactory. 


Conclusion 


Although the experimental results given are 
by no means sutflicient, or absolutely conclusive, 
they indicate that the approach made can, in 
certain cases, be used for the calculation of the 
condensation of binary mixtures with sufficient 
accuracy. 

NOMENCLATURE 

constant in Nusselt’s equation C /(mole /sec)* 
molal concentration moles /em3 
molal specific heat at constant pressure cal mole°C 
specific heat at constant pressure eal /g °C 
diffusion coeflicient em? /see 
acceleration of gravity em /sec* 
= vapour flow moles /sec 
molal enthalpy difference between vapour and 
liquid cal/mole 
= partial heat transfer coefficient cal /°C em? see 
= partial mass transfer coefficient em /see 
flow of condensate moles /sec 
= molecular weight g/mole 
amount of condensate, formed per unit of time 
and area moles /em? sec 

molal flux in r direction - 

molal flux in z direction in 


q 
R 


Re 
r 

S 
Se 
Sh 
T 


1" 


AT 


pP 


Subsi 


= heat flux to the wall cal/em? sec 
== radius of the tube em 
Reynolds number of the vapour no dimension 
= coordinate in radial direction em 
circumference of the tube em 
Schmidt number of the vapour no dimension 
Sherwood number of the vapour no dimension 
== temperature °C 


= temperature of a saturated vapour of a composi- 


tion equal to 7m, C 
= temperature drop from vapour to wall °C 
= mole fraction of more volatile component in the 

liquid phase 
= mole fraction of more volatile component in the 

vapour phase 
co-ordinate in vertical direction (measured from 
the top) em 
length of tube em 
= condensate flow per unit of circumference by 
volume cm*/em sec 
thickness of liquid film em 
related to 
diffusion em 


thickness of “laminar vapour film ” 


= thickness of “laminar vapour film” related to 
heat transfer em 
viscosity g/cm sec 

= heat conductivity cal/°C. em see 


density g/em* 


ripts 
lower boiling component 
= higher boiling component 


' = condensation 


T 
Z 


= Vapour 

= interface 
= liquid 

= conduction 


' = constant or mean composition 


= constant temperature 
= bottom of the tube 
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Abstract 


Under constant drying conditions the j,, factor for heat transport and the jp factor 


for mass transport have been determined for the drying air with four evaporating substances : 


benzene, carbontetrachloride, ethanol and water. 


mentioned substantially j,,; = jp, as wanted. 


Résumé 


It has been found that under the conditions 


Les auteurs ont déterminé, dans des conditions de séchage identiques, les facteurs jj, 


de transfert de chaleur et jp» de transfert massique, relatifs au séchage de lair avee quatre agents 


dévaporation benzene — tétrachlorure de 


carbone 


éthanol et eau. Dans les conditions men- 


tionnées ci-dessus, les expériences montrent légalité entre j,, et jp) comme prévue par la théorie. 


Tue analogy between heat transfer and mass 
transfer, expressed e.g. in the equality of the j 
factor of Cuimuron and CoLpurn [1] for heat 
transfer j, to that of mass transfer jp, for 
geometric equal cases, can be verified in an 
elegant manner in the case of simultaneous heat 
and mass transfer from one body (and therefore 
of exactly the same geometrical form) as occurs in 
drying a moist solid substance in a hot gas. 

In particular, in the constant-rate period of 
drying under constant drying conditions of a 
porous nonshrinking body, the study of the 
validity of the equality mentioned for the drying 
gas is facilitated not only because the surface 
temperature is constant in this period, but also 
by the fact that the resistance to heat flow in 
the gas phase is large, compared with the resist- 
ance in the moist porous body. Therefore, during 
the constant-rate period of drying, not only the 
surface, but the whole body is at a fairly constant 
temperature, viz. the wet-bulb temperature of the 
drying air, because the differences in temperature 
in the body are within the error of determination 
if the bodies chosen are not too thick.* This 
facilitates the measurements, because in such a 
case this temperature can be actually measured 
under drying conditions, and need not be cal- 
culated. Based on this, the present authors have 
undertaken to measure the j,, and jp factor for 


air of four liquids evaporating from a porous 
earthenware body, viz.: benzene, carbontetra- 
chloride, ethanol, and water. 


EXPERIMENTAL 

The principle of the measurement used was that a stream 
of air of known temperature and humidity was passed 
over & moist porous earthenware bloc. The conditions 
chosen were such that the air, in drying, did not change 
appreciably in temperature and moisture content. The 
moist body suspended via a microburette filled with the 
liquid to be evaporated from one of the beams of a balance ; 
the amount of liquid evaporated was then found by 
differences in weight, as a function of the time. If the 
position of the whole setup is rigid enough towards 
horizontal movements, the current of air has no influence 
on the weight measured. 

The setup with the microburette as given was chosen 
to lengthen the constant-rate period. The point of the 
burette was inserted into the porous body and the flow 
of liquid was adjusted in such a manner than it was about 
equal to the amount evaporated. To prevent an eventual 
fall of drops from the body, this was tightly wrapped in a 
rayon gauze. Bare blocs of an originally rectangular 
form of 2 10 em, slightly streamlined at the ends, 
were used. Besides the rate of evaporation, the dry-bulb 
temperatures of the body at several places were measured 
with thermocouples. Moreover, the difference of pressure 
in the measuring section with barometric pressure was 
determined with a differential manometer. 

The air used was conditioned in the following manner 
(see Fig. 1). Moved by a blower, it was first saturated 
with vapour by passing it through a humidifying tower. 
After humidification, the air passed a condenser (cooled 
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* Experimental proof of this fact will be given on another occasion. 
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Fic. 1. Sketch of the drying apparatus. 


with water at constant end-temperature lower than the 
temperature of the humidifying tower), and then a cyclone 
to catch liquid drops. In this way, air of a constant 
humidity (determined at constant speed by the tempera- THERMOCOUPLE 
ture of the condenser) could be obtained. The wanted — 
temperature of the air was then adjusted by passing the 
air through a regulated electrical heater, commanded by 
the entrance temperature of the tunnel. After passing the 
drying tunnel, the air re-entered the humidifying tower. 


During the measurements the humidifying tower was also Yes main SE BA SSSSSY . 
( 
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bypassed. A closed-circuit apparatus had been chosen to 
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prevent great losses of liquid and contamination of the 
atmosphere. If the measurements were carried out with 
the organic liquids, the air was first dried by means of 
silicagel contained in a tower in a bypass, before being 
conditioned. The dew point of the dry air was in all 
these cases 30°F. The metal drying tunnel consisted of a 
calming section and the actual measuring section made 
from copper tube. During conditioning the measuring 
section was bypassed, because it took about 2hr to 
obtain a constant humidity and temperature of the 
air. The amount of air was measured with flow raters 
or a gas meter. The humidity was determined via the 
dew point with an apparatus as sketched in Fig. 2.* 
The long leg of the apparatus was placed in a cooling 

medium, say liquid air, at a temperature far below the Fic. 2. The dew-point apparatus. A = 
dew point to be measured. Therefore heat flows from the — brass cylinder with mirrored surface ; B = 
ebonite isolating piece; C = brass rod ; 


ddlddéddd 


tt antenna 
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* The apparatus has been developed in this laboratory D = metal attaching ring; E = partition 
by C. VAN per Spore. and C, HuyGen on a principle put of case; F = copper swivel; G = ring 
forward by Ir. G. van Tot of the Philips Laboratories at of sponge rubber; H = adjustable air 
Eindhoven. buffer. 
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small cylindrical body A of the apparatus, and as soon 
as its surface has reached a temperature below the dew 
point of the surrounding air, amongst others a con- 
densate will be formed on the mirrored surface of the 
apparatus. A beam of light falling on a photoelectric cell 
via this mirror will then be diffused. 

The interceptance of the beam of light puts in action 
a relay, working on the photocell, by which an electrical 
circuit, in which a small heating element surrounds the 
cylindrical body, is closed. The temperature of the body 
therefore rises and the condensate disappears from the 
The 
described will then repeat with ever smaller and smaller 
The 
mean temperature of the body (measured with a small 


surface, upon which the heating stops. process 


differences in temperature of the cylindrical body. 


thermocouple or with a thermometer) will approach the 
dew point of the surrounding air. The apparatus gives 
dew points with an accuracy of about 0-1 F. However, 
the time to reach equilibrium with sufficient accuracy is 
rather long, (15 minutes) and the instrument therefore can 
only be used when working under stationary conditions, 
As a typical example, in the following Table 1, for one 
experiment chosen at random, the temperatures involved 
are given as a function of the time during stationary 
experimental conditions. 7, and Ty, represent the tem- 
perature at two places in the body. T, is the dry-bulb 


temperature of the air. 


Table 1. 





0 


31-7 
31-6 
17-9 


TC 
TC 
Ty Cc 

Dewpoint 
F 82 





The wet-bulb temperature (T,,) of the air measured in 
a small quantity of liquid in the drving tunnel was 31-7 ¢ 
and therefore equal to the temperature of the body. 


METHOD OF CALCULATION AND RESULTS 


From the amount of liquid evaporated as a 
function of the time the rate of evaporation ¢,, 
can be calculated, The rate of heat transfer dy 
is therefore : $j), = dy, JW, in which IW is the 
enthalpy difference between the vapour at the 
temperature 7’, and the liquid at the temperature 
Tx. The partial coefficient of mass transfer (k) 
and heat transfer (h) can then be calculated with 
the aid of the equations : 


oy = h(T py — Ty) (1) 


Hreermes and W. P. Rincens 


RT 
It has assumed that in 
Nu = ¢(Re)"(Pr)" and Sh 

exponent n }, 
circumstances the air velocity at the interface 


ou (Pw — Pg) (2) 


the correlation 
ec ( Re)" (Sc)" the 
because under the prevailing 


been 


gas body can be assumed to be zero [2].* Therefore 
for j, and j,, the original definitions of Cuiron 
and Co_purN could be used.* 


h ( Pr)? 3 (3) 


Ju (Nu) (Re) (Pr) 1/3 
“p pl 


Jp (Sh) (Rey (Sey 4 


~ (Sc)? 3 ( 1) 


the mixtures 
evaluate equations (3) and (4) have been calcu- 
lated, if need be, 
» T+ 


”» 


The properties of necessary to 


for a mean temperature 


» and for the pressures measured 


and based on numerical data taken from the 


literature. 

For the diffusion coeflicient and the Schmidt 
collected 
Picrorp [3] have been used. 


SHERWOOD and 
For the heat of 
evaporation, the vapour pressures, the density 


number the data of 


and the viscosity, several sources [4] have been 
consulted, 

It has moreover been assumed that the density, 
the specific heat, and the Prandtl number were 
additive for the mixtures. The velocity V has 
been taken as the superficial gas velocity, there- 
fore calculated on the empty tunnel. 

Finally, the attention may be turned to the 
fact that the factor A has not been corrected for 
heat the 
roundings to the bloc, as follows from the equation 


transported by radiation from sur- 


given. Under the conditions chosen this amount 
of heat is small compared with the heat trans- 
ported by conduction and convection. Moreover, 
the calculation of this amount gives an approxi- 
mate figure only. The influence of the radiation 
can be measured by comparing the A factor with 
an analogous factor (/,) for radiation. This 
partial coeflicient of heat transfer by radiation 
can be calculated by means of the equation 

* The present authors are of the opinion that in general 


the j factor should be defined by : j,, = Nu Re Pr” and 
ip = Sh Re" Se, n ranging from 4 to }. 
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ot ‘TW A l 
o( 47,8) | “" 
/ €h “Bw \"x 


')]. (5) 


In this equation, which is valid for small tem- 
perature differences betwene the source and sink 
denotes the arithmetic 


, 
m 


(see e.g. Annoup [5]), 7 


mean temperature of wall and bloc; «, and «, are 


mu 


the emissivities, and A, and A, are the surfaces 


of wall and bloc respectively. The constant o has 
the value: 5-73 10° joules ‘m® see (K)*. 
In the case investigated the ratio A 


u A, ~ 2, 
therefore the emissivity of the wall has a large 
influence on h,, and because the emissivity of the 
red-copper polished wall is small, 4, will conse- 
uently be small. 

m : BIB K, 


0-06, the value of h, 0-7 


Introducing e.g. in equation (5) for T 
for & 0-9; for «, 
Joules m?® see °C, 

When fh, is compared with the values found 
for h (60-27), and the accuracy of the measure- 
the that 


correction of h (about 2%) is not necessary and 


ments seen, conclusion is warranted 
not justified. It may be emphasized that, as has 
been demonstrated by ARNOLD (/.c.) the amount 
of heat transported by radiation, even for small 
temperature differences between source and sink, 
cannot always be neglected, in particular if the 
body to which heat flows is small compared with 
the radiating walls. 

The results obtained are collected in Table 2. 
As can be seen, the variation in the ratio jy jp 


in the various experiments is rather large. How- 


Ju, or J do not seem appropriate, 


ever, no dependence of the Reynolds number can 
be traced, as can be found by plotting jy, /jp 
versus the Reynolds number. It seems therefore 
safe to calculate a mean value for every system. 
These values have been included in the table. 
Including the uncertainty caused by some of 
the data used in the calculation, in particular with 
the diffusion the 
Schmidt number, it can be concluded that, for 


regard to coetlicients and 
the regions and systems investigated, the assump- 


tion jy 
Correlation of the ) factor with the Reynolds 


jp is valid. 


number takes the form 
j c(Rey™~' (6) 
Plotting log 7 against log Re (Re calculated for 
the air in the empty tube), as has been done in 
Fig. 3, enables to calculate ¢ and m by the method 
0-45 and 


It may be observed that in Fig. 3 


of least squares. It appears that ec 
m 0-62. 
some experiments with naphthalene have also 
been included, These experiments have been 
carried out with blocks of pure naphthalene, 
shaped in the same fashion as the earthenware 
bodies. Only the j, factors could be included 
the differences in T, — 7 


small (0-3°C) to enable the calculation of j,, with 


am are too 


because 
enough accuracy. 

Comparison of the data obtained with the 
many found in the literature [3, 6] based on jp, 
It would only 


add to the confusion in this field. The more so, 











j=045R,938 


—+ 








carbon tetrachloride 
carbon tetrachloride 


270 300 320 430 470 300 


—_—-- --—m 10° 


jy benzene 
Ip benzene 
jy «ethanol 

ethanol 


. ip naphthalene 


3. 3. LogJ as a function of log Re. 
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Table 2. 





Water 


Carbontetrachloride 


Benzene 


Ethanol 


12 
3-66 
404 
iS4 
3-07 


414 


mean : 


mean : 


mean : 


mean : 


1-13 
0-96 
0-05 
14 
0-08 
102 
1-02 
144 
0-89 
10 





since the Reynolds number used has not been 


based on the correct parameters. First the super- 


ficial velocity has been used and not a velocity 
with respect to the surface of the block. Secondly, 
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the diameter of the drying section has been taken 
instead of a length parameter from the body. The 
rather peculiar shape of the body prevented the 
authors from doing so with any amount of 
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certainty. Therefore, the only conclusion of some 
importance to be drawn from Fig. 3 is that it 
indicates that equation (6) seems to be fulfilled, 
for the conclusion of which the figure has been 
presented, and that m = 0-62 for the form 
given. Any change in the velocity or the length 
parameter, which will be relatively the same for 
every experiment, does not change this con- 
clusion. 


CONCLUSION 
Under the experimental conditions given, the 


equality of the j,, and the jp factor for the gas 
phase seems fairly well satisfied. 
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NOMENCLATURE 

A = surface m? 
bm = Tate of mass transfer moles /m? sec 
¢y = Tate of heat transfer Joules /m? sec 
T,.» = wet bulb temperature of the air Cc 
Tp = dry bulb temperature of the air *c 
W = enthalpy Joules /mole 
h = partial coefficient of heat transfer Joules /m? °C sec 
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Separation of m-cresol and p-cresol by extractive crystallization 
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Abstract 


The principles of the new separation method, called “ Extractive Crystallization ™ 


are described in detail. This new method has been applied for the separation of cresol mixtures, 
using acetic acid as the extractive crystallization solvent. The results obtained show that this 
hew process can serve as a useful method of separation in the presence of a suitable solvent. 


Résumé 


nommeée “ cristallisation extractive.” 


Les auteurs décrivent en détail les principes d'une nouvelle méthode de séparation 
Ils ont appliqué cette nouvelle méthode a la séparation 


de mélanges de crésols en utilisant l'acide acétique comme solvant d'extraction. 


Les résultats obtenus montrent que ce procédé peut servir utilement pour la séparation en 


présence de solvant approprié. 


INTRODUCTION 


One of the most frequently occurring problems 
in the field of chemical engineering is the separa- 
tion of the components of a mixture. There are 
available several general techniques such as 
fractional distillation, azeotropic or extractive 
distillation, and liquid-liquid extraction, whereby 
such a separation may be brought about. If the 
components of a mixture, however, have com- 
parable volatility and close chemical properties 
so that no selective agent can be found for 
facilitating their separation, none of the above 
methods will be applicable or economic for their 
separation. 

The mixture of m- and p- cresols is such an 
example. Various attempts made in the past to 
separate them from their mixtures by extractive 
or azeotropic distillation have not met with any 
success [14]. These components have been, 
however, resolved by other methods. The com- 
monest method employed is that of treating the 
mixture with a foreign agent which forms 
addition compounds with the individual com- 
ponents. The foreign agent is so chosen that 
either a compound is formed with only one of 
the components or, if the compounds are formed 
with both the components, they may have differ- 
ent physical properties, with the result that they 
can be separated by any of the common methods 
of separation including fractional crystallization. 
For example, BentLey and CatLow [1] used 
benzidine as a foreign agent which formed 


addition compounds with cresols, and these were 
then separated by cooling the mixture to a 
temperature of 95°C (when the p-cresol addition 
compound crystallized out) and subsequent 
filtration. Savirr and Orumer [15] washed the 
benzidine addition compound cake free of 
entrained mother-liquor with an inert solvent 
(benzene), and improved the purity of the 
product from 85°, as obtained by Bentriey and 
Catiow (loc. cit.) to 98%. Among the other 
agents employed the following may be mentioned : 
phenol [10], nitrous acid [8], urea [17], hexamethy- 
lene tetramine [3], 4-picoline [3] and oxalic acid 
[7]. The second method employed is to prepare 
suitable derivatives. which are easily separable 
and then the original compounds won back from 
them. For instance, m- and p-cresols were con- 
verted into their respective sulphonic acids which 
were then separated, taking advantage either of 
their differential solubility in water or of the 
difference in the temperature at which they 
hydrolysed [2]. The cresols were also separated 
by debutylation [19], of the butylated deriva- 
tives. On a laboratory scale, cresol mixtures 
were also separated by liquid-liquid extraction 
using a buffered solution having a pu of 10-8 as 
the extracting solvent [4]. 


An entirely physical method of separation of 
such mixtures into their components is possible 
if the components form eutectic mixtures without 
the formation of solid solutions. For, in absence 
of the formation of solid solutions, one of the 
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components (with only a negligible amount of the 
other component) crystallizes out on cooling till 
the eutectic temperature is reached. The separa- 
tion realizable by such a simple freezing process, 
however, is limited by the composition of the 
eutectic mixtures. However, if the crystallization 
is carried out in presence of a third agent (solvent) 
which forms single eutectic systems with either 
of the feed components, the further separation of 
the components may be possible, and particularly 
so if the ternary eutectic composition is shifted 
in a favourable direction. 

A process of separation of m- and p-cresols 
based on this method (termed hereafter “ ex- 
tractive crystallization ”), using acetic acid as the 
foreign solvent, was tried out, and the present 
paper incorporates the theoretical considerations 
of the above process and the solid-liquid phase 
equilibrium data for the ternary system m-cresol 
p-cresol-acetic acid, 

The term “extractive crystallization ” 
also been used by previous workers to describe 
some of the methods [17] mentioned earlier. 
However, it is felt that this term should exclusively 
be used for a process in which the separation of 


has 
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Fig. 1 (a). 


components themselves is effected by crystalliza- 
tion in presence of a solvent. 


THEORETICAL CONSIDERATIONS 

For discussing the principles underlying the 
process of extractive crystallization, consider a 
ternary system of components A, B, and C, of 
which B and C are the feed components to be 
separated from each other and A is the foreign 
agent. Let this be the simplest possible ternary 
system, viz. the system in which the corresponding 
binaries are all single eutectic systems and in 
which no solid solutions are formed. The nature 
of the solid-liquid equilibrium diagram for such 
a system will be as shown in Figs. 1 (a) and 1 (b). 
Fig. 1(b) consists of a concentration triangle 
ABC with the three binary diagrams laid out 
along its sides. 

In these diagrams, FL (or F’L), KL (or K’L) 
and HL (or H'L) are the binary crystallization 
curves, and they trace out the concentration 
regions AF’LK'A, BH'LF'A and CKLH'C (see 
Fig. 1 (b)) in which the components A, B and C 
will precipitate out in nearly pure forms during 
primary crystallization range. The point L is 
the ternary eutectic point, and its composition 
represents the limit of separation in an extractive 
crystallization process, It is evident from these 
diagrams that if the composition of the overall 
mixture is such that it lies in either of the two 








Fic. 1 (b). 
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triangles K’'LC or H’'LB, both in the primary 
as well as the binary crystallization regions, only 
one of these components and or the solvent will 
crystallize out, and if the crystallization process 
is arrested just before the ternary eutectic point 
is reached, the solid phase will consist of only 
the solvent and one of the feed components, while 
the liquid phase will contain all of the other 
feed component together with the balance amount 
of the solvent and the crystallizing feed com- 


ponent. On the other hand, if the overall com- 


position of the mixture, to start with, is such that 
it lies on either the line CL or BL, right up to the 
eutectic point, only one component (cither C or 
B) will crystallize out, and from the practical 
point of view it would be advantageous to work 


with such mixtures. The diagram also reveals 
that, if the disposition of the ternary eutectic L 
is favourable, the binary eutectic mixtures be- 
tween B and C can be further treated for the 
separation of either B or C in the pure form. For 
example, a required amount of A may be added 
to the mixture H so that a ternary mixture of 
the composition represented by a point P lying 
on the line CL is obtained which on cooling will 
precipitate out the component C right up to the 
ternary eutectic temperature. 

The above considerations also apply for a 
ternary mixture in which one of the three binaries 


E 








is a two-eutectic system accompanied with a 
compound formation as is the case in the m-cresol 
p-cresol-acetic acid system. In such a case, as is 
there will be two ternary 
The point Q lying on the 


indicated in Fig. 2, 
eutectics, R and 8. 
l’A (where VW’ 
binary compound between B and C) may be 


line is the composition of the 
regarded as the eutectic of the pseudobinary 
system between V and A. If it is intended to 
separate B from C, the concentration region 
CRSBC should be avoided. 
minimum points, instead of lying in two different 
triangles (A VB and AVC), may also lie in one 
of them only, as indicated in Fig. 3. In such a 


The two ternary 


case, only one of them is a true ternary eutectic, 
while the other is called a peritectic point. Besides 
these peritectic type of systems, other types of 
systems with solubility gaps in either the liquid 
or solid phase are also possible, but are not 
discussed here as the solid-liquid phase equilibrium 
diagram for the system studied does not belong 
to any of them. 





Cc 


Fic. 3. 


EXPERIMENTAL 
Materials 
Cresols—During the investigation, B.D.H. quality cresols 
(reported as 98-100°, pure) were used after distillation 
in a column packed with glass helices (equivalent to 
twelve theoretical phates). Only the middle cut of the 
fraction distilling over in the temperature range of 
201-202°C was used, which usually had the following 
constant melting points for the two cresols : 
11-50°C, and p-cresol 34:5°C. 
Acetic acid—The Merck quality acetic acid was first 
refluxed for about two hours with a calculated amount 
of acetic anhydride, and thereafter the middle cut of the 


m-cresol 
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fraction distilling over at 118°C was collected and had a 
melting point of 165°C. 


Determination of binary solid liquid equilibrium diagrams 


The apparatus for the measurement of disappearance 
points consisted of a freezing point tube of about 1} in. 
in diameter which carried a mechanical stirrer and an 
accurately calibrated thermometer. About 25g of the 
mixture of a known composition was charged into it, and 
thereafter it was surrounded by an unsilvered vacuum 
flask (about 4 in. i.d.) which carried an alcohol-solid carbon 
dioxide freezing mixture. The distance between the bottom 
of the tube and the surface of the mixture determined 
the rate of cooling. Where supercooling with the mixture 
was excessive, it was necessary to freeze the whole mixture 
and then warm it slowly by taking it away from the 
freezing mixture, this being effected by the lowering of 
the Dewar flask. The warming rate was adjusted to about 
0-5-1-0°C a minute to enable the accurate determination 
of disappearance points. The whole mixture was illumi- 
nated by an electric bulb outside the Dewar flask to 
enable careful observation of either the disappearance or 
appearance points. (Appearance points could only be 
determined with acetic-acid-rich mixtures.) 
mixture, the process was repeated three or four times and 


For each 


the mean value of the readings was taken as the correct 
temperature at which the phase change occurred. The 
solid-liquid phase equilibrium data so obtained for the 
binaries m-cresol—p-cresol, m-cresol-acetic, and p-cresol 


acetic acid are plotted in Figs. 5 and 6. 


Determination of binary eutectic curves 


To locate the binary crystallization curves in the ternary 
system, m-cresol-p-cresol-acetice acid, the disappearance 
or appearance points were determined for a large number 
of synthetic mixtures over the entire concentration range 
using the same apparatus as described above. Ternary 
mixtures investigated were usually prepared by adding 
an accurately weighed quantity of the pure component 
into a mixture of the other two components of known 
composition, and thus covering the mixture compositions 
lying in a vertical section passing through one of the 
corners of the temperature composition prism. Starting 
with different binary mixtures, the whole concentration 
range was thus investigated. The disappearance points of 
the ternary mixtures lying in any one vertical section were 
plotted against their respective compositions, expressed 
on the basis of the percentage amount of the third com- 
ponent in the original binary mixture. The minimum 


point or points in these curves represent the points on 
the binary crystallization curve or curves which are cut 
across by the vertical sections. 

These minimum points when plotted on an equilateral 
concentration triangle (Fig. 7) trace out the binary 


crystallization curves. These curves tend to converge 
into the two points representing the compositions of the 
two ternary eutectic points. By working with mixtures 


near the possible intersection point of the binary eutectic 
curves the eutectic ternary composition was determined. 
The composition so determined was verified by noting 
whether it froze or melted at a constant temperature like 
a pure component. 


Separation of the two phases in equilibrium 


The apparatus used for separating the two equilibrium 
phases is indicated in Fig. 4. It consisted of a jacketted 
equilibrium tube E (about 1} in. in diameter) in the top 
part of which was fused a quickfit socket joint in which 
fitted a brass cone carrying a stirrer assembly. The lower 
part of the tube carried a sinter glass plate within it and 
a bottom connection through which the liquid phase was 
withdrawn. The desired temperature in the equilibrium 
tube was maintained by circulating cold aleohol through 
the jacket round it. The cooling and the temperature 
control systems were arranged as shown in the figure. The 
seepage of the liquid through the sinter disc while it was 
being maintained at the equilibrium temperature was 
prevented by filling the free volume below the dise with 
mercury, the level of which could be raised or lowered with 
the help of the mercury reservoir W as shown. When the 
equilibrium conditions were reached, the mercury was 
withdrawn through the side arm and the mixture was 
filtered to the receivers placed below by applying vacuum 
to them. First the liquid phase was withdrawn and 
collected in one of the receivers B. The crystals remaining 
behind (together with the entrained liquid) was then 
allowed to melt by stopping the alcohol circulation and 
was thereafter withdrawn into the other receiver. 

This apparatus was found to be very satisfactory from 
the points of view of quick filtration, ease of operation, 
and minimum transfer losses. The filtration was normally 
over within a fraction of a minute except when the solid 
to liquid ratio was high and the liquid phase very viscous, 
these conditions being usually encountered at temperatures 
lower than — 15°C. Even for such cases, the filtration 
time was never higher than 2-3 min. 

In some of the runs, attempts were also made to remove 
as much as possible of the entrained mother-liquor from 
the crystals by displacement washing. Most of the common 
organic solvents were found miscible with the cresols and 
hence the sodium chloride brine saturated at the equilibrium 
temperature was employed for washing. This was, how- 
ever, not tried for all the runs as it introduced fresh 
complications such as the preferential washing of acetic 
acid by the brine and an increase in the time cycle of 
filtration, especially at the low temperatures. 


Analysis of binary and ternary mixtures 


For the estimation of acetic acid, an accurately weighed 
amount of the ternary mixture (about 1 g) was diluted to 
100 ml. An aliquot of 25 ml was then titrated with N/10 
sodium hydroxide, using as an indicator a mixture of 
equal amounts of 0-1°, solutions of bromthymol blue and 
cresol red, which gave a sharp colour change from yellow 
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R= Alcohol circulation 
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Fig. 4. 


to dark green at pu 7:5. The method was found capable 
of giving results accurate to within + 0-1°,. 

For the estimation of cresols, about 25 ml of the same 
ternary mixture were distilled in a fractional distillation 
assembly until the residue in the still as tested was com- 
pletely free of acetic acid. The fraction coming over 
thereafter contained cresols, and although some cresols 
were lost during the removal of acetic acid, it was found 
that their relative proportions remained unchanged. This 
is because the vapour pressures of these two cresols and 
their mixtures are found to be identical over the entire 
concentration range {12}. 
found to be ineffective in preferentially lowering the 
partial vapour pressure of any one cresol in the mixtures 
[18]. Disappearance points of the acetic-acid-free mixture 
were then determined with the apparatus previously 
described, and its composition was read from the binary 
temperature-composition diagram as determined. The 
accuracy of this method, checked with synthetic mixtures, 
was found to be within 2°,. 


Presence of acetic acid also is 


Discussion AND RESULTS 
A. Binary solid-liquid equilibrium diagrams 


These diagrams for the three binary systems are 
represented in Figs. 5 and 6. Both m-cresol-acetic 
acid and p-cresol-acetic acid are single eutectic 
systems, while the binary m-cresol-—p-cresol is a 
two-eutectic system, accompanied with a com- 
pound formation. It may be observed that the 
phase behaviour of both the cresols with acetic 
acid is very similar, particularly in the acetic 
acid rich region wherein the freezing temperatures 
of various binary mixtures are almost identical 
and determined by the amount of cresol (irres- 
pective of which one) in the original mixture. 
The m-cresol-p-cresol system has been studied 
previously by Dawson and Mountrorp [6] and 
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Present iwestigations hence only the disappearance points which were 
Data by Dawson found to be always reproducible were relied on. 
and Mountford The method of thermal analysis also was not found 
Data by Stevens : , ee ee » Aime : " 
ond tiituin satisfac tory in view of the difficulty of super 
cooling of the cresols. 

The melting points of the pure components 
employed and the eutectic data as determined, 
are presented in Table 1 given below. 











B. The ternary system 


The concentration triangle for the ternary system 
as determined is represented in Fig. 7. In this 
diagram R and S are the two ternary eutectics, 
while the curves N’R, K'R, F’S and J’S are the 
four binary eutectic crystallization curves. QR 
and QS can also be considered as pseudo-binary 








crystallization curves. These curves were located 
by plotting the composition of the minimum 
point or points obtained in the primary crystal- 





lization curves for various vertical sections 
judiciously chosen to cover the entire concen- 
tration range. Along the curves N’R and F’S, 
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Fic. 5. System m-cresol—p-cresol. 








also by Stevens and Nicke.s [20] whose is a 
more recent work. The melting points of the 
cresols used by us and of the cresols used by the 
latter were almost identical and, in spite of this, 
a fair divergence is observed in their and our 
results. This divergence can therefore be explained 
only on the basis of the method employed for the 
determination of freezing points. In view of the 
difficulty of supercooling with cresol mixtures, 
in the present work only the disappearance points 
were determined with about 25 g of the sample 
which was kept mechanically stirred. As against a System p-cresol 
this, Stevens and Nicke.s used 5-10g of the =. —{ acetic acid 
sample, employed seeding and hand-stirring with o System m-cresol 
a wire loop, and considered the maximum t— acetic acid = 
temperature reached during freezing as the correct ; | 
temperature. Such a procedure for a mixture 10 20 30 40 50 60 70 80 90 100 
and one which supercools considerably even when *JeAcetic acid ——+ 

seeded was not found satisfactory by us, and Fic. 6. 





+ dun 
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acetic acid with either p-cresol or m-cresol will 
co-precipitate, while along the curves A’'R and 
J’'S binary compound V’ with either p-cresol or 
m-cresol will crystallize out. Similarly, along the 
pseudo-binary crystallization curves QR and QS 
acetic acid will co-precipitate with the compound 
Vv’. Thus, if the purpose is to separate the two 
cresols, the concentration region CRS BC should 
be avoided. The composition of m-cresol in the 
commercial mixtures is usually of the order of 
60°, which lies between the composition (on the 
acetic-acid-free basis) of the two ternary eutectics, 
and therefore it is evident that just adding 
acetic acid alone to it and crystallizing will not 
effect the separation of the two cresols. This 
will be possible only when the ternary system 
with the solvent is a peritectic type in which 
both the true eutectic and the peritectic points 
lie on either left or right of the line 4P where the 
point P represents the composition of the com- 
mercial mixture. 
of either of the two cresols which can be recovered 


Also the percentage amount 


in the pure form from the mixture P will depend 
upon the distance of the nearest eutectic point 
from the line AP. The greater this distance, the 
greater will be the temperature and consequently 
the concentration range over which the separation 
can take place. Similarly, it is obvious that the 
separation of one of the components in the pure 


K 


of AV no Ms 
0 0 20 38 40 SO 60 7 8 90 WO 
p-Cresol m-Cresol 





Fig. 7. System m-cresol-p-cresol-acetic acid. 

form from a binary eutectic will take place only 
if the composition of the ternary eutectic (expressd 
on the solvent-free basis) is different from the 
composition of the original binary eutectic 
mixture ; and the larger will be such a separation, 
the greater the difference in the ratio of the 
amounts of feed component in the binary and 
the ternary eutectic mixtures. As will be seen 
from the concentration triangle, the two ternary 
eutectics in this case do not lie very far away 


Table 1. 





Melting point, C 


Composition, ©, m-cresol 
First Eutectic 
Melting point, “C 


Composition, °,, m-cresol 
Binary compound 
Melting point, “¢ 


» 


Composition, ©, m-cresol 


Second Eutectic 


Melting point, °C 


Melting point, °C 


Present 
investigation 


Dawson and 
Mountrorp [6] 


STEVENS and 
NIcKELs [20] 
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from the lines K'A and J'A, and as such not 
much separation of the two binary eutectics 
(K' and J’) in the presence of acetic acid is to 
Attempts were made to resolve 
two 


be expected. 
the 
cresols, but the separation theoretically expected 


above two eutectic mixtures for the 
being very low and the experimental errors 
involved being comparatively large, no conclusive 
results were obtained. 

However, a large number of runs were carried 
out on the separation of mixtures lying outside 
the 
results obtained are recorded in Table 2. On the 


basis of the equilibrium temperature-concentra- 


these eutectic compositions, and some of 


tion triangle as determined, and the lever rule, 
the per cent recovery expected from each run at 
any given temperature has been calculated, and 
these figures are also included in the same table 
for comparison. To explain the procedure followed 
in calculating the per cent recovery of one of the 
components in its pure form at any given tem- 
perature of crystallization, consider the mixture 
of the composition 16°, m-cresol, 65-6°,, p-cresol 


and 18-4°, acetic acid (Run No. 4, Table 2) 


maintained at an equilibrium temperature of 


15°C. By interpolation, it is found from the 
disappearance point data of mixtures that the 
composition of the liquid in equilibrium with the 
15 € 17-5°, 


p-cresol, and 21.5°, acid. 


solid at should be m-cresol, 
610°, The 
composition of the equilibrium liquid phase as 
collected analysed 17.3°, 61-5°, p- 
cresol, and 21-4°, acetic acid, indicating fairly 


acetic 
m-cresol, 
good agreement between the calculated and the 


observed values, particularly in view of the fact 
that the composition of the equilibrium liquid 


phase was arrived at by interpolation of the 
These 


disappearance-point data on mixtures. 
data for various runs carried out are included in 
column 3 of Table 2. In column No. 4 are given 
the compositions of the residue (as recovered after 
filtration) consisting of the solid phase and the 
entrained liquid. If it is assumed that the liquid 
entrained has the same composition as that of 
the equilibrium liquid phase, the composition of 
the solid phase as recovered will be obtained by 
drawing a straight line through the points 
representing the compositions respectively of the 


residue, the equilibrium liquid, and the feed 
mixture. Following this procedure, the com- 
position of the solid obtained is indicated in 
Fig. 8, and the line drawn through the points 
representing the compositions of the feed mixture 
and the equilibrium liquid (as calculated) illus- 
trates the divergence between the two results, 
and it will be found that it is well within the 
experimental error. On the basis of the com- 
positions of the solid phase arrived at as above, 
the the 
ponent, calculated on the basis of its amount in 


percentage amounts of desired com- 


© Observed values Acetic acid 


x Calculated values 


the feed mixture, are also calculated, and are 
represented in column 6 of the same table. These 
figures for the mixture taken for example are 
21.0 and 20.5 respectively. Divergence in some 
other cases is rather high and is mostly to be 
attributed to the error involved in the estimation 
of the correct composition of the equilibrium 
liquid, and to a smaller extent to the inaccuracy 


involved in the analysis. 


CONCLUSION 


The solvent employed acetic acid — to separate 


m-cresol and p-cresol from their mixtures by 


extractive crystallization, is not particularly 


suitable, as the composition of the two ternary 


eutectics on the solvent-free basis are not 
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very different from those of the corresponding extractive crystallization process is a potential 


binary eutectics between the two cresols. This method of separation if a suitable solvent. is 


work, however, does conclusively prove that the found. 
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Book reviews 


High 
Vol. I. 
F. DD. 


1056. p. 


speed aerodynamics and jet propulsion, 
Thermodynamics and Physics of Matter. 
(Editor), Oxford University 
100s. (in U.K.). 


Rossini Press, 


S12 xviii. 
Tuts work, to be completed in twelve volumes (of which 
the first is here reviewed), is intended to “ bring together, 
in convenient, highly accurate, and well organized form 
the basic information that has evolved in the last decade 
in the related fields of gas dynamics, aerodynamics, com- 
bustion and jet propulsion.” This first volume in the series, 
edited by Rossint, brings together “ the basic principles 
of thermodynamics and physies of matter, presented from 
the standpoint of their applicability to a better under- 
standing of high speed aerodynamics and jet propulsion.” 


Rossini, covers the fundamentals of 
In section B, Herzrecp and GrirrinG 
discuss quantum mechanics, including molecular structure, 
Hirscu- 


statistical 


Section <A, by 
thermodynamics. 
bond energies, and activation energies ; while 
FELDER, Binp and Sporz 
mechanics and the kinetic theory of gases. Section C, by 


CURTISS, discuss 
Breartie, gives a complete treatment of the thermodynamic 
properties of real gases and their mixtures. In section D, 
Birp, and Sporz 
transport properties of gases. In section FE, Rice covers 
Section F, by and 
Brinkiey, presents the theories of liquids and liquid 


HinscureLper, CURTISS, discuss the 


critical phenomena. RICHARDSON 
In section G, Ewaxp discusses solids and solid 
Section H, by Herzrenp, covers relaxation 
In section I, EsrerMANN presents 
a discussion of gases at low pressures, including molecular 


solutions. 
solutions. 


phenomena in gases. 


beams. In section J, Curtiss discusses the thermodynamics 
of irreversible processes. 


The task of presenting such a wide range of subjects 
in so short a compass is no easy one, but the reviewer is 
not convinced that the method used here, of compiling 
contributions from acknowledged leaders in their fields, is 
the best, or indeed a good, solution. The defects are those 
that derive from the method of presentation : there is 
unevenness of treatment, and there is repetition. Both 
may be illustrated by comparing the treatments of wave 
mechanics in sections B and G. There are sections that 
are simply condensed or rewritten from parts of other 
books, and it is not always clear what has been gained 
in the condensation. There is a lack of unity which is 
perhaps inevitable. 

To comment on sections in detail is probably to reflect 
principally the reviewer's interests. Section A is an 
adequate introduction, although the ideas of temperature 
and of scales of temperature have been better expounded 
elsewhere. Perhaps the most valuable section will be found 
to be that (section C) on the thermodynamic properties 


of real gases, where Beatrice collects a large number of 
equations for the evaluation of actual properties in real 
cases. It may prove to be unfortunate that Rice, in his 
otherwise excellent discussion of critical phenomena, has 
done so on the basis of the assumption of a critical region, 
an assumption that is novel and is still debated. The 
article on relaxation phenomena in gases is succinct and 
well documented. That on the thermodynamics of 
irreversible processes (written in 1951) is now out of date, 
and fails to point out the potential importance of this 
field. 


A feature of the book 
bibliographies to each section, 


is the lists of references and 
although the 
considers that more attempt should have been made to 
bring these up to date. There are notable absences, for 
example in the field of thermodynamics, which reflect the 
long interval between the writing of the section and the 
date of publication. A defect is the absence of reference 
tables showing the symbols used and their meanings ; this 
detracts from the usefulness of the book as the work of 
reference it is intended to be. 

The standard of production is high, and so is the price. 


reviewer 


H. P. Hurcutson 


TREYBAL. 
$9.50 


Roserr E. 
> 666 pp. 


Mass-Transfer Operations. 
McGraw-Hill Book Company, 1955. ix 
(£3 11s. 6d.) 


Tue author has written an integrated account of the 
chemical engineering applications of mass transfer, and 
the work is so arranged that the essential similarities 
bet ween involving the fundamental 
mechanisms are stressed. The particular features of the 
book which commend it are the chapters dealing with 
diffusion in solids and with adsorption. The last-mentioned 
which also treats ion exchange as a special application of 
adsorption, gives a useful introduction to a subject which 
is not included in most general textbooks. The title of the 
book is, however, rather misleading, since this is not an 
advanced book on mass transfer, Apart from the items 
mentioned above, no material is included which one would 
not expect to find in a reasonably advanced textbook of 
chemical engineering, and some of the operations are not 
considered in quite the detail which might have been 
expected in a book specializing in mass-transfer opera- 
tions. Thus, only passing reference is made to the penetra- 
tion theory of Hicsre, and nowhere is a method of solution 
given for the unsteady mass-transfer equation. Air 
humidification is quite well covered, though the method 
of calculation of the condition of the gas phase is not 
given. In the absorption chapter, no reference is made to 


processes same 





Book reviews 


simultaneous absorption and chemical reaction, and the 


newer centrifugal absorbers are not mentioned. Again, 


multicomponent distillation is not treated, and the book 
therefore falls short of the requirements of the honours 
undergraduate student in this respect. Liquid-liquid 
extraction is covered very fully, as one might expect in 
view of the author's interest in the subject, and the 
chapter on drying is comprehensive, though one would 
little 


moisture movement in the capillary 


mechanism of 
The des- 


criptive section of the chapter on leaching is good, but 


have liked to see a more on the 


Spaces, 


a little more could profitably have been said on the 


methods of calculation for systems involving counter- 


current washing. One of the most useful features of the 


book is the large number of worked examples, and the 


inclusion of a wide selection of problems at the end of 


each chapter for the student to work through himself. 
Altogether, the book is a well-written exposition of the 
subject and gives a somewhat novel approach to several 
of the operations discussed, Except in the field of adsorp- 


tion, however, it cannot be said to add a great deal to 
the scope of existing textbooks. 
M.W.P. 


Journal of the Imperial College Chemical Engineer- 
ing Society, Volume 8 Ltd., 
Redhill, 1954. 112 pp. 12s. 6d. 


Love and Malcomson, 


InN the eighth volume of the journal, the high standards 
set by the preceding volumes have been maintained. This 
publication is not simply a house journal of the Imperial 
College, but contains a number of original contributions 
of considerable scientific importance by past and present 
members of the chemical engineering department. In 
some cases the work referred to is still in progress, and 
the journal thus serves to give advance information on 
subjects to be reported more fully at a later date. Among 
the subjects included are freeze drying, rheology, measure- 
ment of gas solubility at high pressures, and silicones. 


PR, 





Announcement 


Victoria, Australia, offer £1,000 Fellowship in Chemical 
Engineering 


Tue Fellowship worth £1,000 (Australian) is offered in 
Chemical Engineering by the British Memorial Fund, 
established in Victoria, Australia, by public subscription. 
The Fellowship is a gesture of “ loyalty, gratitude and 
affection to the British people in recognition of their 
“ role in saving civilization in the Second World War.” 


It provides for ten months’ study in Victoria for the 


selected candidate and covers travelling and living expenses 


during the period. 


Accommodation will be arranged by the Fund. Wives 
and children are welcome but no extra expenses can be 
paid to cover them. The upper age limit for both men and 
women is 35 vears, and applicants must be of British 
stock of at least ten years’ residence in the United King- 


dom. 


The Fellowship is open to workers in any aspect of 
chemical-engineering research or  chemical-engineering 
education who could benefit by their participation in the 
development of this science in Australia. The following 
topics are listed as a guide to the types of investigation 


work which might be undertaken by the Research Fellow : 


1. Economic and production problems of chem- 
ical manufacture in Australia. 
Process efliciency — methods of investigation. 
Gas turbines in chemical industry. 
Water purification and conservation in arid 
areas. 
Chemical engineering in extensive agricul- 
ture. 


6. Dynamics of road surfaces subject to heavy 
traffic. 
Free radical mechanisms in chemical indus- 
trv. 
Furnace design and combustion efficiency. 
Heat and mass transfer to and from the 
atmosphere. 
Zone refining with particular reference to 
solutions of inorganic salts. 
Dimensional factors in heterogeneous cata- 
lytic process. 

The necessary arrangements and general direction of 
research work will be the responsibility of 1.C.1. Australia 
and New Zealand Ltd., who have their headquarters and 
major research facilities in the State of Victoria. 

Simultaneously with the Fellowship in Chemical 
Engineering, a Fellowship in Natural History and two 
others where applicants may choose their course of study 
or research are being offered by the Fund. 

Members of the London Selection Committee are : 


Col. the Hon. W. W. Lecoarr — Chairman 
Lady ALBEMARLE Lord HUNTINGFIELD 
Lord BAILuret Sir Vincent Trewson 

Sir ANGus GILLAN 

Dr. J. F. Foster 
Association of Universities of the British 

Commonwealth). 


(Secretary 


Full details can be obtained from Col. the Hon. W. W. 
Leecarr, Agent-General for Victoria and Chairman of the 
British Memorial Fund, London Section Committee, 
Victoria House, Melbourne Place, Strand, London W.C.2. 
Closing date for the receipt of applications — 15 September 
1956. 





